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ABSTRACT 


Nap-of-the-earth (NOE) flight in a conventional helicopter is ex- 
tremely taxing for two pilots under visual conditions. Developing a 
single pilot all-weather NOE capability will require a fully automatic NOE 
navigation and flight control capability for which innovative guidance and 
control concepts have been examined. Constrained time— optimality provides 
a validated criterion for automatically controlled NOE maneuvers if the 
pilot is to have confidence in the automated maneuvering technique; this 
is one focus of this study. A second focus has been to organize the stor- 
age and real-time updating of NOE terrain profiles and obstacles in 
course-oriented coordinates indexed to the mission flight plan. A method 
is presented for using pre-flight geodetic parameter identification to 
establish guidance commands for planned flight profiles and alternates. A 
method is then suggested for interpolating this guidance command informa- 
tion with the aid of forward— and side-looking sensors within the 
resolution of the stored data base, enriching the data content with real- 
time information, and combining the stored and sensed data for real-time 
display, guidance, and control purposes. A third focus of this study 
defined a class of automatic anticipative guidance algorithms and neces- 
sary data preview requirements to follow the vertical, lateral, and 
longitudinal guidance commands dictated by the updated flight profiles and 
to address the effects of processing delays in digital guidance and con- 
trol system candidates. The results of this three-fold research effort 
offer promising alternatives designed to gain pilot acceptance for automa- 
tic guidance and control of rotorcraft in NOE operations. 
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SECTION I 


BACKGROUND AND INTRODUCTION 


A. BACKGROUND 

The U.S. Army and helicopter manufacturers are vitally interested in 
automatic nap-of-the-earth (NOE) flight. The more perceptive of these 
believe that most of the guidance and control problems associated with 
Light Helicopter Experimental (LHX) and AH-64X rotorcraft operations have 
at least interim solutions using a mixture of manual and automatic re- 
sources. There is some promise that the range of automatic operations can 
be extended further to include automatic NOE flight. Before this can be 
accomplished, however, there are at least three major issues that must be 
resolved: 


1. Interpolating within the resolution of a stored NOE data base 
with the aid of forward— and side-looking sensors, enriching 
the data content with real-time information, and combining 
the stored and sensed data for real-time display, guidance, 
and control purposes. 

2. Development of path and attitude command signals from this 
data array that are appropriate to command safe NOE flight of 
a three-dimensional helicopter (in contrast to a point mass 
helicopter) . 

3. Pilot acceptance of automatic NOE flight. 


In the past, the path toward full automation of dangerous mission- 
critical elements that could, in principle, be accomplished by the pilot 
has been difficult and fraught with subtle misunderstandings between the 
piloting and engineering communities. Automatic landing is perhaps the 
best example of this. Ultimately successful systems required the develop- 
ment not only of guidance and control systems that could do the job with 
acceptable performance and fail-operational reliability, but they also 
required : 
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® Displays presenting status and back-up command information 
with which the pilot could monitor the progress of opera- 
tions, the automatic equipment, etc., with the option of 
instantaneous takeover 

• Control laws that, to a very large extent, mimicked the op- 
erations of a pilot when following the same guidance 
information. 

The last point above is indeed quite subtle, because it implies that a key 
factor in pilot acceptance of automatic equipment to fly complex maneuvers 
is that "it flies the airplane as a pilot would." 

For the automatic equipment to fly the helicopter as a pilot would, an 
automatic system is required that emulates not only the error-correcting 
features of ordinary feedback control- processes but also a system that has 
predictive and anticipative properties. In terms of the classical modes 
of pilot control behavior, the automatic equipment should incorporate so- 
called "pursuit" and, for some maneuvers, "precogniti ve" elements in 
addition to the "compensatory" (error correcting) feedback control system 
behavioral aspects common to both automatic and manual control systems. 
These higher-order modes of control behavior are not easily mechanized in 
typical guidance systems, because they require prediction and 
anticipation. The pilot develops the higher-level behavioral modes by 
discerning more cues (e.g., via preview) than are conventionally available 
to monometric single-quantity sensors and by virtue of operations as a 
dual-mode controller when the input command is appropriate. Thus the 
pilot in higher level control is very sensitive to the total command as 
well as any errors, while most automatic equipment is largely error 
driven. 

B. INTRODUCTION 

Contemporary technology has provided automated guidance systems for 
up-and-away flight between stored waypoints using inertial, Doppler, 
Loran, TACAN, VOR/DME, VORTAC, Omega, and Global Positioning Systems 
(GPS). Contemporary technology has also provided independent automatic 
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terrain following/terrain avoidance (TF/TA) and course-keeping systems for 
low altitude flight using forward looking radar. 

Piloted guidance technology for rotorcraft, furthermore, has now pro- 
vided passive FLIR and SLIR sensors and image converters. This technology 
could be adapted with additional processing of the converted images, e.g., 
as suggested in Ref. 1, in order to develop an innovative form of automa- 
tic guidance and control which is compatible with piloted guidance. We 
shall, however, adopt a somewhat different approach here, based on the 
proposal in Ref. 2. 

Automatic nap-of-the-earth (NOE) operations based on stored knowledge 
of terrain, obstacles, threats, and targets will require advancements in 
two related technologies. The first is the creation of accurate geodetic 
bases for navigation, threats, and targets in any theater of operation. 
The second is the development of storage and display media for real-time 
use of the data bases in flight. The geodetic bases must cover wide areas 
in order to allow for versatile mission planning; this will require data 
storage media more compact than now available (Ref. 3). 

For flight planning and in-flight uses, however, the pilot usually 
prefers that the terrain profiles, obstacles, threats, and targets be 
presented in course-oriented coordinates for monitoring automatic guidance 
as well as for manual backup guidance. Although the adoption of course- 
oriented coordinates will inherently contribute significant compression of 
terrain data storage requirements for a specific mission, additional data 
compression is desirable based on the properties of the terrain itself. 
This issue is addressed in Refs. 4 and 5 and in Section II, "Task I: 
Storage and Updating of Terrain Profiles and Obstacles," which describes 
our innovative method for organizing the storage and realtime retrieval of 
terrain profiles, sensed obstacles, threats, and targets in course- 
oriented coordinates indexed to a defined flight profile. 

Some of the consequential guidance and control data processing issues 
associated with the automation of following defined flight profiles in NOE 
operations are discussed in Section III, "Task II: Automatic Guidance for 
Following Flight Profiles," and in Section IV, "Task III: Automatic 
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Guidance for Aggressive NOE Maneuvers* 1 ’ Both of these sections have 

supporting appendices* 

Section V presents a summary of pilot-centered considerations for 
monitoring automatically controlled NOE flight* These considerations 
include automatic-manual control response compatibility, status displays, 
and decoupled control responses, all of which are essential for pilot 
acceptance of automatic NOE flight. 

Finally, Section VI presents a summary of the principal findings for 
each of the tasks that constitute this phase of the study. 
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SECTION II 


TASK Is STORAGE AND UPDATING OF TERRAIN PROFILES AND OBSTACLES 


On future U.S. Army helicopters, such as the Light Helicopter- 
Experimental (LHX) , the entire scout or attack mission may be accomplished 
in instrument meteorological conditions (IMC) with a single pilot. Con- 
sequently, if the LHX concept of a single-place attack (or even a scout) 
mission is to be viable, many of the elements in a typical sequence would 
have to be handled by automatic equipment. 

Virtually all functions of scout/attack rotorcraft depend on maintain- 
ing the planned flight profile as a function of time. Other mission 
functions — such as threat warning and countermeasures, target detection 
and acquisition, and weapon delivery and fire control — are intimately 
related to the flight guidance function. Therefore, it is logical to 
consider automatic navigation and flight control as "core functions" for 
automating the flight guidance functions. 

In an automated guidance system for rotorcraft, a computer-assisted 
briefing system should augment or supplant the verbal transmission of 
mission data. An important part of an automated guidance system would be 
a prepared portable storage medium that would permit automatic loading of 
all mission— required data into the rotorcraft^s guidance system. The 
mission data required would include (from Ref. 3): 

• Routes, waypoints, initial points, targets 

® Terrain profiles 

• Obstacle types and locations 

• Threat types and locations 

• Target characteristics 

• Weapons complements. 
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It would be necessary to update automatically this pre-flight data base 
with information from scouts and other airborne warning and control sys- 
tems via a secure data transfer link* 

Automatic NOE operations based on stored knowledge of terrain, ob- 
stacles, threats, and targets will require advancements in three related 
technologies. The first of these is the creation of accurate geodetic 
bases with sufficient resolution for navigation, threats, and targets in 
any theater of operation. The second technology is the development of 
storage and display media for real-time use of the data bases in flight. 
The geodetic bases must cover wide areas in order to allow for versatile 
mission planning; this will require data storage media more compact than 
are now available (Ref. 3). The third technology is the interpolation 
within and enrichment in real time of the stored data base with a sensed 
data base of obstacles and threats acquired from forward- and side-looking 
sensors. 

For flight planning and in-flight uses, however, the pilot usually 
prefers that the terrain profiles, obstacles, threats, and targets be 
presented in course-oriented coordinates for monitoring automatic 
guidance. (Airborne radar, lidar, and infrared sensors can also be 
directed to scan in course-oriented coordinates.) Part 1 of Table 1 shows 
that U.S. Army standards for scout/attack terrain flight training require 
navigation within 1500 m left or right of the preselected course. Thus 
one innovation which we have begun to investigate is to organize the 
storage and real-time retrieval of terrain profiles, sensed obstacles, 
threats, and targets in course-oriented coordinates which are indexed to 
the defined flight profile. Besides condensing the data storage 
requirements, this arrangement will simplify the coordinate conversion and 
identification required for compatible processing of real-time data from 
on-board forward-looking and side-looking terrain/obstacle/threat 
sensors. This processed real-time data will be applied subsequently in a 
coordinate estimation algorithm to update the defined flight profile and 
the guidance commands derived therefrom for use by the automatic guidance 
algorithms. The automatic guidance algorithms will employ estimates of 
predicted vertical, lateral, and longitudinal deviation commands required 
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to follow the updated path, course, and schedule, respectively, with 
feedforward command compensation. Preview of both stored and sensed 
terrain and obstacles is therefore essential. Novel approaches to setting 
both compensatory and pursuit feedforward preview distances are 
recommended subsequently for minimizing the guidance error in the presence 
of atmospheric disturbances which may be correlated with the terrain and 
obstacles. 

A. A METHOD FOR PREFLIGHT GEODETIC PARAMETER 
IDENTIFICATION TO ESTABLISH GUIDANCE 
COMMANDS FOR PLANNED FLIGHT PROFILES 


This method involves the assumption of a mathematical structural model 
of flight profile guidance commands with undetermined coefficients, and it 
determines the coefficients of the required guidance command strategy by 
recursively matching the model to the sensed obstacle and terrain data 
using a running least-squares estimation technique (i.e., a multiple lin- 
ear regression). Initially, the coefficients of the model are defined by 
the guidance commands for the a priori flight plan. (In flight, the 
guidance commands will be updated in real time from a sensed data base 
using a procedure to be discussed in a subsequent subtopic.) 

The guidance command vector, R(jX , t ) , will be a function of the state 
vector, X_. Both R_ and X_ will have along-course, across-course, and eleva- 
tion components. The guidance command vector, R_(X_,t) will also be a 
function of events in real-time coordinate, t, detected by the terrain and 
obstacle sensors or communicated by data link. In the mathematical model 
that we have investigated, R(X_>t) will be represented by a combination of 
polynomial and harmonic functions of along-course and across-course co- 
ordinates, i.e., a vector of truncated Taylor and Fourier series. 

The numerical analysis concept central to the procedure consists of 
applying a running least-squares estimation technique using a moving data 
"window." A specific example of least-squares estimation occurs in curve- 
fitting problems where it is desired to obtain a functional form of some 
chosen order that best fits a given set of discrete measurements. Thus 
the least-squares problem that we have investigated consists of 
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correlating a dependent variable, z, e.g., the elevation component of 
R(X,t), with one or more elements of an independent variable vector, X. 
It is assumed that z and _X are related by an expression of the form 

z = f lCl + f 2 c 2 + ... + f j Cj ... + f m c m (1) 

z = [f 1§ f 2 , ... ] [cj , c 2 , ... cJ T 

or z = F_c_ (2) 

where f^ , f 2 , ... are selected variables from the state vector X or ex- 

plicit functions thereof, and the Cj's are constant coefficients relating 
y and f j • In Eq. 2 is a row— vector of the f j , and c is a column— vector 

of the c j . For example, z_ = F_ c_ will include the following forms for 
Fourier and Taylor series in one dimension. 

Fourier Series 


z(A,C) = c 


m 


0 


or 


+ Sc. cos (\. - E ) where the c. are real 

1=1 J J 1 J 


z( A) 


E c exp ( i A) where the c. are complex 
j =-n 3 J J 


Taylor Series 


z(x) = 


P 

d^ + Ed, 
o . . k 
k=l 


^ x ~ x ref ) 
k! 


where d. 


x<k) ( x ref) 


For a terrain model, however, z is also a function of across-course 
coordinate y, i.e., z = z(x,y^) as shown in Table 2 for truncated Fourier 
and Taylor series representations. 
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TABLE 2. TERMS FROM TRUNCATED FOURIER AND TAYLOR SERIES 

FOR TERRAIN MODEL 


FOURIER SERIES 

Assume the Fourier series in the form of a sum of cosines and sines 
with undetermined coefficients and given periods for f kffl : 


Define: 


C km^km 


2c. 

~L 


km 

2nnrx 

- - - .L. 

i 

T 

c 

JS 

rg 

2mirx 
./ 

jm * 0 

— C08 
X 

™ 1 T 

L 

X 

OAU 

L 

X 

L 

X 

|k - 1 


1, 2, 3, 


.. M 
K 


Height of 


Across 



- z(x, y k ) - c k0 + f kl c kl + f k2 c k2 + f k3 c k3 + c km 

Estimate the kth set 

of c km U8in 8 y " yk 

and 2M values of x 
to determine f(x,y k ) 

Model requires 2M*K values of z corresponding to 2M values of x at K 
values of y to estimate c km (y k ) and d km (y k ) 

TAYLOR SERIES 


In particular situations it will be necessary also to incorporate 
terms from truncated Taylor's series: 


f (k+1 )(m+l ) C (k+l )(m+l ) 
f (k+2)(nH-2) C (k+2)(m+2) 


^ x_x o^ C (k+l ) (m+1 ) for linear slope 
( x-x 0 ) 2 ° (k+2 ) (m+2 ) for curvature 


10 



An example of the array of undetermined coefficients for coordinate 
zCx.y^) for a terrain model (and for a corresponding guidance command 
model) is shown in Table 3, with m harmonics along-course and k coordin- 
ates across-course. Examples of the specific number (1 + 2m)k of 
numerical coefficients required to represent a 3 km x 30 km terrain model 
of a route having a particular x-y coordinate resolution are listed in 
Table 4. The present study investigated the representation of profiles 

from a 5 km x 5 km digitized terrain model furnished by NASA and having 
10 m resolution shown in Fig. 1. 

The number of coefficients of the terrain model (and corresponding 
guidance command model) required for a truncated Fourier series represen- 
tation can be estimated in advance by preprocessing sections of the 
terrain for the intended route with a finite Fourier transformation 
(FFT). The significant harmonic content for guidance can be estimated 

from the effective spatial bandwidth of the terrain over the route to be 
followed by the rotorcraft. Figure 2 provides a convenient graph for 
relating spatial bandwidth of the terrain or route to the temporal 
bandwidth of a guidance and control system at three forward speeds. At 
20 kts, spatial bandwidth in radians per decameter (rad/dm) and temporal 
bandwidth in radians per second (rad/sec) are virtually identical 
numbers. The sections of terrain in Fig. 1 for which FFTs were obtained 
exhibited effective spatial bandwidths of about 0.2 rad/dm. Thus we can 
expect a typical guidance and control system to follow an offset profile 
over these samples with good precision at 20 to 40 kts forward speed. 

The STI proprietary computer program for obtaining FFTs from digitized 
samples of data is called FREq uency Domain Analysis (FREDA). The NASA-STI 
computer program for parameter identification using a running least 
squares estimation is called Non-Jjitrusive Parameter _Identif ication 
program (NIPIP). These acronyms will recur in subsequent text. Before we 
present some of the results using FREDA and NIPIP, we shall digress to 
summarize the steps in the procedure used in NIPIP to obtain least— squares 
estimates for the undetermined coefficients in the terrain model. 
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TABLE 3. TERRAIN MODEL ARRAY OF c km (y k ), d fcm (y k ) 
AS SHOWN WITH (1+2M)K COEFFICIENTS 
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TABLE 4. EXAMPLE OF ROUTE 
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Procedure for Obtaining Least-Squares Estimate 
For Undetermined Terrain Model Coefficients 

If there exist sets of discrete measurements for z and F from the 
terrain and obstacle sensors and navigation systems and the difference 
between the actual guidance profile required to avoid terrain and ob- 
stacles and the output of the assumed model given by Eq. 1 is given by v, 
then 


z = F c + v 
n — n — n 


where £ is the estimated value of £ given by Eq. 1, and the subscript n is 

a 

used to denote each set of y and F_. It is assumed that there exists a c 

which will relate all values of z„ to the F . 

n — n 

To arrive at an estimate for c^ at least m sets of measurements are 
required (where m is the number of degrees of freedom in Eq. 1). The 

A 

solution for c_ is found by adjoining these sets of discrete measurements 
as follows. 


“ — 
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Z 1 

Z 2 
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v i 

V 2 
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II 

• 

£ + 

m 

• 

Z N 




. 

• 

V N 

_ _ 


— _ 


_ 


(3) 


where N is the number of measurements and N > m. Equation 3 can be writ- 
ten more concisely in a matrix notation as 


z = H c + v 


(4) 


where _z, H, and _v are defined by inspection. Note that H is an N by m 
matrix. 
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The derivation of the least-squares estimate for the unknown 
coefficients, c_, is obtained by minimizing the sum of the squares of v 
where 



(5) 


The complete derivation is given in Ref. 2. 

A 

The recursive least-squares estimate for c is given as 


N 

£ - ( 2 

n=l 


T 

F F 
— n — r 


^-1 


N 

E 

n=l 


T 

F z 
— n n 


( 6 ) 


where N is the number of data points up to the given point in time where 
the sum is truncated. 

The following features of this solution can be noted. First, the main 

computational task consists of updating summations of products of the m by 

m matrix F_ F_. The only storage requirement is that these summations alone 

be saved, not all of the accumulated data as with a batch processing 

scheme. Second, the most complex computational task is the inversion of 
N 

the E F F matrix. 

. — n — n 
n=l 

This completes our digression to describe the way in which the least- 
squares estimate is performed by NIPIP (more details on NIPIP can be found 
in Refs. 6, 7, and 8). We shall now return to describe the results of 
using FREDA and NIPIP to model terrain using Fourier series descriptors. 

B. RESULTS OF TERRAIN MODELING USING 
FOURIER SERIES DESCRIPTORS 

Recall that the purpose of this effort is to evaluate various ways of 
arriving at a terrain model (and corresponding offset guidance command 
model) that uses a Fourier-series-like structure, but with a reduced 
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number of components. The example terrain data base consists of a set of 
altitudes over an x-y grid. Figure 1 (on p. 14) shows a contour plot for 
the example data base. This data set is available on a MICRO-VAX II as 
discrete altitude in units of decameters (dm) stored as an array having a 
discrete x-y grid of 512 x 512 dm. Six slices along y - 110 , 185, 199, 
200, 201, and 440 dm were unpacked and stored as IBM PC files for subse- 
quent analysis using FREDA and NIPIP. STI's JFrequency CONv ersion program 
FC0N was used to produce the binary input file format for entry into 
FREDA. 

FREDA can provide results only at frequencies that are harmonics of a 
cycle per run length. NIPIP is more general, in that it can use non- 
harmonics in its fitting process, albeit at a cost of larger matrices and 
slower execution in storing all of the cross-products required. Thus, our 
procedure is to use FREDA to determine bandwidth and provide insight into 
those spatial frequencies having predominant power in order to guide the 
selection of a reduced number of frequency components for use in NIPIP. 

I. FREDA Results 

Figure 3 shows plots of the spatial history and power spectral density 
(PSD) of elevation along the slice at y * 200 dm. The spatial history 
includes the average value that was removed during the PSD analysis. The 
logarithmically-spaced bin selection technique was chosen so that the 
first 24 PSD results (Fig. 3) would have a bandwidth of one cycle per run 
length (expressed in rad/dm) and center frequencies at multiples of the 
first one, which is 2 tt/ 512 = 0.0123 rad/dm. Above the first 24 
components, the PSD is binned using an approximate log-spacing technique 
up to the maximum available frequency (tt rad/dm in Fig. 3). Figures 4 and 
5 show results for y = 440 and 110 dm, respectively. While the detailed 
spatial histories are different, there is a second-order roll-off in power 
(first-order roll-off in amplitude) above the effective bandwidth of about 
0.2 rad/dm in all three PSDs. 
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Mean Terrain Height 150.192 dm 
Root Mean Squared Height 54.3 dm 
Mean Squared Height 2946.9 (dm)^ 
Effective Spatial Bandwidth 0.2 rad/dm 


Figure 3. FREDA Output for y ■ 200 dm 
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Figure 4. FREDA Output for y - 440 dm 
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2 at Y = 110 decanetere 



Mean Terrain Height 119.009 dm 
Root Mean Squared Height 29.6 dm 
Mean Squared Height 875.47 (dm)^ 
Effective Spatial Bandwidth 0.2 rad/dm 


Figure 5. FREDA Output for y ■ 110 dm 
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2. NIPIP Results 


The NIPIP program was adapted to the IBM PC in a standard version. A 
special version was then set up to allow it to identify the amplitudes of 
a sum of sines and cosines to represent terrain data. Some additional 
enhancements were needed for dynamic range scaling of the matrices to 

prevent overflows that can occur for a long data window and a large number 
of sines and cosines. These enhancements reflect the results of an 
insightful review of the similarity of the NIPIP and FREDA formulations 
for the special case of sine waves that have an integer number of cycles 
in a data window or run length. For this case, the determinant of a key 
NIPIP matrix is bounded by N*(0.5*N*G*G)**NSC, where N is the number of 

data points, G is the amplitude factor, and NSC is the total number of 

sines and cosines. In the future, this relationship can be incorporated 
into NIPIP to scale the calculations automatically. 

Figure 6 shows the output file from NIPIP with 20 harmonics, and it 

compares some of the peak amplitudes of the sine (AS^) and cosine (AC^) 
components with those found using FREDA. The results are virtually iden- 
tical. Figure 7a compares the 20-harmonic fit to the y = 200 dm slice 
(the smoother curve is the fit). The data compression ratios are shown in 
each caption. 

In Fig. 7a, the resulting extreme deviations of the approximation are 
within 10 dm with three exceptions: the initial valley at X = 0 is 

under-estimated by 12 dm, the valley at x = 225 dm is under-estimated by 

15 dm, and the mesa at x = 512 dm is under-estimated by 20 dm. Obviously, 

under-estimation of the peaks is unsafe, but the extreme deviations of 
peaks can be more confidently predicted by sensed information; whereas, 

the valleys are usually shielded from forward-looking sensors by interven- 
ing peaks. Thus the procedure is at least on the safe side (with respect 
to terrain) in under-estimating the valleys, because the affected flight 
profile commands, based on the stored knowledge of terrain, will be up- 
dated in real time by coordinate identification using sensed obstacle and 
terrain data. This will provide the capability to interpolate within 
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a) First 20 Harmonics. Data Compression Ratio 512: [1+20(3)] , i.e., 
512:61 or 8.39:1 



b) First Five Harmonics Plus Two Non— Harmonics at 0.1 and 0.14 rad/dm. 
Data Compression Ratio 512: [1+7(3)], i.e., 512:22 or 23.3:1 



c) First Five Harmonics. Data Compression Ratio 512: [1+5(3)] , i.e., 
512:16 or 32:1 


Figure 7. Effect of Number of Harmonics on Terrain Fits for y ■ 200 dm 





geodetic data bases having a coarse quantization of terrain, such as the 
Defense Mapping Agency (DMA) 10 dm data base. 

Figure 7c shows the effect of reducing the frequencies to the first 
five harmonics; whereas, Fig. 7b shows the effect of adding two non- 
harmonic frequencies (0.1 and 0.14 rad/dm) above the first five. These 
preliminary efforts show that about 20+ frequencies are needed if the 
peaks are to be accurately described. Finally, Fig. 8 shows the fit for 
the first 206 dm using the same seven frequencies used in Fig. 7b. This 
case is included here to give some initial insight into what can be 
achieved with non-harmonic frequencies selected from inspection of the 
PSD, since the fit is clearly superior for the same region. 

These results suggest that the techniques examined can achieve practi- 
cal data compression ratios between 8:1 and 10:1 for the samples of 
vertical profiles from the terrain data base in Fig. 1. The resulting 
recommendations for compressing flight profile data storage requirements 
are summarized in Table 5. 

This completes our discussion of the results of terrain modeling for 
the purpose of representing off-set guidance commands using Fourier series 
descriptors. In the next subtopic, we shall describe a procedure for 
updating and modifying the stored flight profile guidance commands in real 
time using a data base derived from forward- and sideward-looking sensors 
and having a (navigation) reference system in common with that for the 
stored data base. 

C. REAL-TIME MODIFICATION OF THE FLIGHT PLAN 

TO ACCOMMODATE SENSED TERRAIN AND OBSTACLES 

Figure 9 illustrates a procedural flow diagram for updating and modi- 
fying the stored guidance command flight profiles in real time to provide 
unexpected obstacle avoidance. Two data bases having a common naviga- 
tional reference system (which provides state vector 20 are identified at 
the top of the diagram: the stored flight profiles (vector R J1 ) at the 
upper left and the sensed profile of terrain, obstacles, and threats, 
together with offset bias requirements for safety at the upper right 
(vector R ). 
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Figure 8. Seven Frequency Fit to the First 206 dm Portion of the 
y - 200 dm Slice with Non-Harmonically Related Frequencies. 
Data Compression Ratio 206: [ 1+7(3) ] , i.e., 206:22 or 9.36:1 
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TABLE 5 


RECOMMENDATIONS FOR COMPRESSING FLIGHT PROFILE 
DATA STORAGE REQUIREMENTS 


Preprocess flight profile with FFT program, e.g., FREq uency J)omain 
Analysis (FREDA) Program 


Identify mean bias, slope, and curvature in planned flight profile 


Identify special patterned features such as step functions, square 
waves, trapezoidal waves, triangular waves, and sawtooth waves, which 
require extraordinary spatial frequency bandwidth to represent 
harmonically 


Estimate predominant spatial harmonic content in planned flight 
profile after identifying, defining, and removing mean bias, slope, 
curvature, and special patterned features which can otherwise be 
represented by truncated polynomials and simpler unique functions 


Define criteria for adjusting Am to smooth raw spectra 


Select only those frequencies among u> m (m • 1, 2, . .., m) that 

represent predominant power in the FFT of the flight profile 


Omit many of the higher frequencies with low power, and use non— 
harmonically related frequencies in selected regions 


Use only the selected frequencies in the vector of Fourier series with 
undetermined coefficients 


Reduce the number of undetermined coefficients to be stored 


Modify parameter identification program (NIPIP) for automatic self- 
scaling 


Consider shorter run lengths, although this will require transition 
logic to blend the ends during reconstruction 
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Profile 
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Command and Re- 
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Flight Profile 


Dash 


Select 
Dash/ Stop 
Command 


Figure 9. Procedural Flow Diagram for Unexpected Obstacle 

Avoidance in Real Time 
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In principle, it is necessary to compare R with in real time and 
resolve any conflicts with the planned flight profile by automatically 
selecting the indicated lateral or vertical offset maneuvers and returning 
to the planned flight profile where possible. If lateral or vertical 
offset maneuvers will not remove a conflict between R and R , it may be 
necessary to accelerate or decelerate or stop and let the pilot select 
another flight plan. The automatic obstacle avoidance maneuvers will be 
selected from the constrained time-optimal repertory to be discussed sub- 
sequently in connection with Task III. A practical method for comparing 
R^ and R^ is described briefly in the following paragraphs. 

A three-dimensional safety margin envelope is defined mathematically 
in the form of a cylinder which encompasses the extremities of the rotor- 
craft with room to spare, governed by safety margins. This cylindrical 
envelope is mathematically centered on the planned course and vertically 
offset profile. When the rotorcraft is hovering, the cylindrical envelope 
is centered on the present position of -the rotorcraft. When the rotor- 
craft is translating, an hemicylindrical envelope bisected by the 
direction of the inertial velocity (Vj) advances mathematically ahead of 
the rotorcraft's present position a distance Vjx p , where x p is the preview 
time interval of the anticipative trajectory coupler discussed in 
Section III with respect to Task II. 

If the sensed profile of terrain, obstacles, and threats (i.e., the 
sensed data base) mathematically penetrates the anticipative hemicylinder, 

a conflict exists along the planned flight course-and-prof ile. This 
conflict must be detected in real time by coordinate identification and 
comparison over the leading surface of the anticipative hemicylinder using 
sensed terrain and obstacle data. It remains to develop in Phase II of 
this study a rationale and procedure for correcting the conflict and 
returning to the flight plan safely by means of transition logic between 
the anticipative trajectory coupler (discussed in Section III) and the 
constrained time-optimal maneuvers for avoiding the obstacles (discussed 
in Section IV). 
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SECTION III 


TASK II: AUTOMATIC GUIDANCE FOR FOLLOWING FLIGHT PROFILES 


A. GUIDANCE AND CONTROL DATA PROCESSING 

Automatic processing of stored and sensed terrain and obstacle data 

will involve sampling, interpolation, extrapolation, smoothing, and com- 
parison of spatial and temporal series of data coupled with priority logic 
for use of the processed data by the automatic guidance algorithms. An 
example of a first level of priority logic for unexpected obstacle avoid- 
ance using discrete maneuvers was illustrated in the flow diagram, Fig. 9 
(p. 28). The automatic guidance algorithms will employ estimates of pre- 
dicted (i.e., previewed or time-advanced) vertical, lateral, and 

longitudinal deviation commands required to follow the desired path, 
course, and schedule, respectively, with feedforward command compensa- 
tion. Preview of both stored and sensed terrain and obstacle data is 

therefore essential in order to compensate for the inevitable (and pre- 
dictable) processing delays and rotorcraft response lags. 

A novel approach to setting the (stored and sensed terrain data) pre- 
view distance has been developed in which the best command-following 

flight guidance and control system is synthesized, and its equivalent 
closed-loop time delay, t e , for low-frequency inputs is determined. The 
data preview distance, Dpe , is then set to read (stored and sensed) ter- 
rain at Dp e = Vt e ahead of the rotorcraft in the direction of travel at 

velocity, V. (Alternatively, in the time domain, the data preview inter- 
val T = D /V = t . ) 
p pe e 

Precision NOE guidance and control, however, imposes a number of other 
requirements on data processing in real time. Consider, for example, the 
vector block diagram representing a multiloop guidance and control system 
in Fig. 10 which is subject to response-correlated disturbances, 
8D( s ) / 3M( s ) . Each vector and matrix will be represented in operational 
rotation as a function of the Laplace transformation operator, s. (The 
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various parts of the block diagram that have been addressed by tasks in 
Phase I are identified in Fig. 10 with Arabic task numbers rather than 
Roman.) The stored flight plan (SFP) and compressed data for the guidance 
command R(s) based on the terrain model investigated in Task 1 is shown at 
the extreme left-side of the figure. The command R(s) is applied (with 
preview) to the pursuit compensation (Task 2) at the top of the figure. A 
filtered command Rf(s) is applied to the compensatory error-correcting 
feedback loop including error processing and equalizing element Yg and 
measurement filter with injected measurement noise N(s) added to 

response M(s). A constrained time-optimal program (Task 3) for rapid 
response (RR) maneuvers if error E is greater than threshold E T is vested 
in element Y R ^. The summed control output C(s) from the three levels of 
guidance, viz., compensatory, pursuit, and programmed rapid response, is 
applied to the controlled element at the right side of the diagram. The 
preferred controlled element for this application will be a decoupled 
velocity command-position hold flight control system-and-rotorcraf t . 

When each vector consists of only one component, the system error 
relationships reduce to the following conventional transfer functions 


e(s) = 


[w r (i - 


3d-, 

3m' 1 


W Y Y lr(s) - fl 
m c p J _ 


- Y 


3d 


)W n(s) - W Y d(s) 


A(s) 


and the system response relationships become 


m(s) 



W Y Y n(s) + Y d(s) 
m c e d 


A(s) 


Stabilization and bandwidth requirements vested in A(s) are subject to 
alteration by response-correlated disturbance gradients 3d/ 3m. 
Characteristic singularities of A(s) are given by the roots of 

A(s) = {l + W Y Y - Y, = 0 
1 m c e d 3m J 
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Processing delays T g (computer frame time, data skewness, algorithmic 
delays) in Y g and higher frequency lags represented by effective delay t c 
in controlled element Y c can be at least partially compensated for their 
effects on stability by lead compensation (or prediction) in W m and wholly 
compensated for their effects on response delay by command preview in 
W r . [Approximate the measurement prediction filter by the ideal form: 
W m = exp(x m s) and the command preview filter by the form: W r = exp(x r s)] 

Response-correlated disturbance gradients 3d/ 3m act to increase or 
decrease the effective open-loop gain of the guidance and control system 
and may even be so great as to compromise stability in the presence of 
delay or to compromise bandwidth regardless of delay. For example, if 



Y = K expf-x s) ; 
e e e ' ’ 



the closed-loop guidance bandwidth and stability will be governed by the 
characteristic roots of 


A(s) = s(s + -j-— ) + K K exp [ (r -t -t )s] - A = 0 
v T ' ce r LV m c e ' J 3m 


Examples of the complex root loci without and with delay appear in 
Figs* 11a and lib, respectively* 


That part of the guidance error (e) associated with commands 
(rf = W r r) may be reduced by the ideal pursuit guidance adjustment 



A(s ) = s(s + 1/T C ) + K c K e exp[(x m - t c - T g )sj- A(3d/3m) - 0 


A) Without Dela> 


Closed-loop roots 
if A(3d/3m) < 0 


co - /K K - A(3d/3m) 
n c e 


Closed loop design with 
no wind shear 


to - /K K 
n c e 


£io - 1/2T 
n c 


Closed-loop roots 
if A( 3d/ 3m) > 0 


B) With Dela\ 


-1/2T 


Closed-loop roots if A(3d/3m) < 0- 


Note low 
margin of 
stability 


Closed-loop design with no ■ 
wind shear 


Closed-loop roots if A(3d/3m) > 0 


Figure 11. Effect of Response-Correlated Disturbance Gradients (e.g., 
Wind Shear) on Closed-Loop Roots for Height or Position Guidance 
and Control Without and With Processing Delay 
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If 3d/ 3in = 0 or if 3d/ 3m is estimated by 3d/ 3r t Y by Y , 

— c c 

and Y^ by Y^ such an adjustment is quite practical for 

K 

c 


C c 

which represent the physical translation dynamics of the attitude 
stabilized, velocity-command-augmented, and decoupled rotorcraft. 
[Denominator factor (s+l/T c ) in Y c represents (augmented) heave damping 
(z-axis), sway damping (y-axis), surge damping (x-axis), or yaw damping.] 

The practical adjustment requires . a pursuit guidance command 
Cp(s) - Ypr(s), which becomes 



f-T s ) ; and Y , 
^ c ; c 


s ( s + ?-) 


c (s) = 
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1_ 
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t 8 ( 


s + 


-) - A M] gS r(s, 
c m 


W 


This pursuit adjustment requires a weighted linear combination of the 
acceleration, velocity, and displacement of the guidance course— and— 
profile command r(s). Commanded velocity must, however, be limited in 
accordance with propulsion constraints (e.g., rotor torque and speed 
limits, blade-stall limits, best-climb rationships between collective and 
cyclic blade pitch), attitude and heading limits near terrain and 
obstacles (e.g., main and tail rotor strikes), and control authority. 
Commanded acceleration must also be limited in accordance with rotor 
flapping limitations, attitude limits, and control rate limits. 
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The human operator will adopt at least the first two terms, 
[sr(s)/K c T c ] and [s r(s)/K c ], of this pursuit adjustment in following a 
fluctuating curved course if sufficient preview of the course slope and 
curvature are available to the operator (Ref. 9). (The operator infers 
the first derivative of the course command from the visible slope and the 
second derivative, from the visible curvature.) In addition, the operator 
may adopt the third term, ( A/K c ) ( 9d/ 3r) , if he has knowledge of or prior 
experience with the wind shear field on the particular course being 
flown. If incorporated in automatic guidance, this practical pursuit 
adjustment will be consistent with piloting technique, given sufficient 
visual preview and will therefore enhance pilot acceptance of automatic 
guidance in following NOE flight profiles with precision. A stored data 
base for the flight profile will provide a practical basis for this 
pursuit command in automatic guidance. 

If this practical adjustment is incorporated, the transfer function 
for the part of the closed-loop system response to command r(s) becomes 


Pursuit 

Adjustment 



Compensatory 

Adjustment 



Pursuit Adjustment for 
/ Wind Shear 


m(s) 

r(s) 


fa /W Y + 
L r me 


W Y - [W Y ,/W Y ) 3d/ 3r lY 
re ^ r d m c J J < 


1 + 


W Y Y - 
m c e 


Y, 3d/ 3m 


The effect of the response-correlated disturbance gradient 3d/3ra can be 
approximately compensated by introducing estimated command-correlated 
disturbance gradient 3d/ 3r = 3d/ 3m in the pursuit adjustment, and the 
overall command response fidelity and delay can be compensated by 
introducing preview in command filter W r = exp(x r s), where T r = x m , if 
W m = exp(T m s). A stored data base for the flight profile will provide a 
practical basis for this preview requirement, where the pursuit data 
preview distance will be VjT r » V^x in the direction of travel at 
inertial velocity V... 
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If the pursuit adjustment is not incorporated or fails, the transfer 
function for the part of the closed-loop system response to r(s) reduces 
to that for compensatory guidance: 


m(s) _ 
r(s) 


W Y Y 
r e c 


1 + W Y Y - Y. 

m c e d 3m 


The compensatory response delay in Y e Y c , 


i.e. , x e + t c , can 
approximately compensated by introducing a different preview 
command filter W r = exp(x r s), where x r = T e + t c . 


now be 
in the 


A significant degree of immunity from wind shears and turbulence is 
provided by automatic velocity— command— position— hold guidance and control 
systems (Ref. 10) whose low order effective controlled element forms are 
given in Table 6. With the exception of the forms for the hover turn, the 
forms for the other maneuvers are valid at translational velocities 
typical of NOE operations. The revised form for heading regulation when 
the translational velocity is not zero is given at the end of Table 6. 
Each of these controlled elements has the generic form 


Y 

c 



(7) 


for a displacement response to a velocity command. If we postulate this 
generic form for Y c , representing vehicle displacement response to a 
velocity command, the ideal pursuit feedforward guidance function 
Y p = [ exp( T m s) ] /Y c which operates on guidance command r, becomes 


Y = f exp ft + 
p L ^ m 



( 8 ) 


[exp(x ra + x R )s] 
K 

c 


{s 2 + (2 ecu - Us + [ 
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PRACTICAL EXAMPLES OF LOW-ORDER EFFECTIVE CONTROLLED ELEMENTS FOR 
FIVE NOE MANEUVERS IN ROTORCRAFT 
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This Ideal 


T m + T k of 
velocity 


guidance adjustment requires prediction over time interval 

2 

weighted commanded acceleration s r/K c , weighted commanded 



K 

c 


sr 


and weighted washed-out commanded displacement 






If the damping ratio £ = 1, the weighted washed-out commanded displacement 
becomes 



O) 

n 


-T~) 


sr 


K 


(. +f) 

K 


and the weighted commanded velocity becomes 


(2u) 

^ n 



K 

c 


sr 


If, in addition to £ = 1, w n * 1/T K best design practice regardless 

of damping ratio), the controlled element in Eq. 7 reduces to 


Y 

c 


K e 
c 


■V 



(9) 


for displacement response to a velocity command. The weighted washed-out 
commanded displacement Is unnecessary (except for improving immunity to 
wind shear), and the weighted commanded velocity reduces to sr/l^T^. The 
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required weighted linear combination of predicted commanded acceleration, 
velocity (and washed— out displacement, if needed) can be readily derived 
with a low noise level from the stored flight profile format discussed 
previously as part of Task I (Refs. 4 and 5), since continuous functions 
with continuous first and second derivatives are used to represent the 
required flight profile(s), whereas, the discrete maneuver logic in Fig. 9 
will be invoked to cope with unexpected obstacle avoidance. 

B. FIELD OF COVERAGE FOR SENSED TERRAIN 

AND OBSTACLES 

Other characteristics of the stored and sensed data which affect pro- 
cessing and must be defined are field of coverage, resolution or 
quantization, and update rate. Relationships among horizontal field of 

coverage, preview interval, velocity, and bank angle in both coordinated 
level turns and uncoordinated level side stepping have been examined in 
11. Two relationships are plotted in Fig. 12 for a fixed preview 
interval of 3 sec. The upper figure is for coordinated level turning; the 
lower, for level side stepping. Both relationships are practically the 
same for half-fields of coverage less than 30 deg. A 60-deg half-field 
angle (120 deg field of coverage), on the other hand, will accommodate 
banking at 60 deg for a 2 g level turn at 50 kt and a 2 g level side step 
at 30 kt. Reference 11, however, suggests that 60-deg bank angles are 
rarely used in NOE operations. Instead, reduced bank angles are likely at 
lower speeds — possibly in accord with the "hypothesized NOE boundary" in 
the lower graph in Fig. 11 corresponding to the 60-deg half-field angle 
down to a bank angle of 20 deg. The cut-off at a bank angle of 20 deg 
represents a typical hovering bank angle requirement based on needs for 
wind-proofing and has nothing to do with horizontal field of data 
coverage. Larger bank angles than those implied by the "hypothesized NOE 
boundary" are probably not used, because the resulting horizontal 
accelerations cannot be effectively managed in close quarters. Thus, if 
the preview interval is 3 sec, it is likely that a horizontal field of 
data coverage greater than 120 deg can seldom be used in NOE operations 
involving forward flight. (Clearly, if automatic side stepping maneuvers 
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a e* 










from hovering flight are required in exceptional instances, the defined 

half-field angle of coverage must exceed 90 deg by virtue of a side- 

looking sensor.) 

Conversion of the horizontal half-field angle of coverage to the half- 

width distance, W, of terrain coverage gives 2W = Tig tan . If the 

' P S 

preview interval is 3 sec and the bank angle does not exceed 60 deg, the 
required half-width of terrain coverage to initiate the sidestepping man- 
euver in forward flight will be 76 m (251 ft) left or right of the 

preselected course at commencement of the maneuver. This half-width, W, 
is only about 5 percent of the Army standard 1500 m, listed in Part 1 of 
Table 1 for navigating left or right of the preselected course in terrain 
flight training. Thus there is ample margin in this example for contin- 
uing the maneuver or increasing the preview interval while retaining the 
60 deg bank angle limitation. 

C. EFFECT OF GUIDANCE AND CONTROL DELAY 

ON HEIGHT CONTROL WHILE TRAVELING IN 

"DOLPHIN" MANEUVERS 

The first and most fundamental effect of discrete data processing to 
be considered in a digital guidance and control system is time delay, the 
cause of which may be vested in a particular algorithm, in the overall 
computational frame time or update rate, in data acquisition by a scanning 
and sampling sensor, or in data conversion from digital to analog form and 
vice versa. 

Another form of delay is peculiar to data multiplexing. When data is 
acquired sequentially, converted from analog to digital form and vice 
versa, or when discrete data is passed between computers serially, some 
data will inevitably be "skewed" (i.e., delayed) in time with respect to 
other data. Thus there will be a progressive "staleness" associated with 
skewed samples of data; sometimes the effect of this staleness can be 
significant, although it is usually much less than the processing delay 
discussed above. 
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A third form of delay is called computation "frame slip" and is as- 
sociated with data transfer between two or more independent asynchronous 
digital processors. Each independent processor may have a slightly dif- 
ferent frame time from the others in the network. For example, one 

processor might have a frame time of 40 ms and another, 39 ms. Every so 
often (about every 1.56 sec), the faster processor will have executed one 
more full frame than the slower processor. This constitutes the effect 
called "frame slip," which can produce undesirable transport delay jump 
phenomena in the dynamic system. 

Effects of data sampling, quantization, and control roughness or in- 
tersample ripple in digital guidance and flight control systems can be 
treated by means of an innovative direct digital design procedure evolved 
and applied in Refs. 12 through 20. 


The effects of delays involving processing and higher frequency 
rotorcraft lags on the closed-loop bandwidth of height control while 
traveling in "dolphin" maneuvers will be first illustrated using numerical 
characteristics of the UH-1H helicopter at 20 kts, sea level, 8000 lbs 
weight, with a mid-c.g. location. Subsequently, it will be possible to 
normalize the time delay, the required open-loop compensation, and the 
closed-loop bandwidth of height regulation in terms of the heave damping 
and collective control effectiveness of the rotorcraft and thus to provide 
more general design requirements to counteract the effects of time delay 
in any similar rotorcraft height, heading, or position guidance and 
control system. 


Figure 13 presents a complex root locus showing the effects of varying 

the height-to-collective control high-frequency gain with a constant 

lead compensation 1 /T t = 1.13 rad/sec and a constant value of the 

L h 

airspeed-to-pitch attitude control high frequency gain 

K u = -0.01 rad 9 c /ft/sec designed to provide a closed-loop speed regula- 
tion subsidence 1 /Tq^ = 0.3 rad/sec. There is no delay represented in the 
results of Fig. 12. Closed-loop characteristics are shown for two values 
of K^, 0.0245 and 0.0923 in/ft/sec. The closed-loop damping ratio in both 
cases is at least 0.7; the larger closed-loop undamped natural frequency 
is 1 rad/sec, and the smaller is about half of the larger value. 
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Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant A M (s) of Height and Airspeed Control Equations 

While Traveling 


A”(s) = s(s + 1/Tg )(s + 1/T Q ) - Z 6 K^(s + 1/T ) 

1 2 c h 


• < 8 - \ + l x « >-h K + - * + l x S >-h ) Z J1 

c c c c 

' Is + 1/T*) 


+ K (X - g)s[s - Z + (X Z /X - g)l 
u v a L w v w a a 


’—a 1/T* -f 

7 A 0 


- U + 1/T 0 ) 


Figure 13. Effects of Varying Height-to-Collective Control Loop Gain 
with Lead Compensation 1/T, = 1.13 rad/sec and Airspeed-to-Pitch 

Attitude Control Gain K u = -0.01 rad 0 c /ft/sec 



The effect of a 0.5 sec time delay involving processing and higher- 
frequency lags is shown in Fig. 14. The lead compensation remains at 
1.13 rad/sec, and the airspeed-to-pitch attitude control loop gain remains 
the same at -0.01 rad 9 c /ft/sec. At the lower value of 

» 0.0245 in/ft/sec, the closed-loop characteristics remain approxi- 
mately the same. The damping ratio is 0.69, and the closed-loop undamped 
natural frequency is 0.54 rad/sec. At the higher value of gain required 
to achieve an undamped natural frequency of 1 rad/sec, the damping ratio 
is reduced to 0.54 with 0.5 sec time delay. Figure 15 shows that the 
closed-loop damping ratio can be restored to about 0.7 with an 0.5 sec 
delay by reducing the frequency of the lead compensation, 1/T^ > to a 
value of 0.89 rad/sec and reducing the high frequency gain to a value 
of 0.074 in/ft/sec. 


The effect of a 1.0 sec time delay involving processing and higher- 
frequency lags is shown in Fig. 16. The lead compensation is back at 
1.13 rad/sec, and the airspeed-to-pitch attitude control loop gain remains 
the same at -0.01 rad 0 c /ft/sec. If = 0.0205 in/ft/sec, the closed- 
loop damping ratio is 0.61, and the closed-loop undamped natural frequency 
is 0.50 rad/sec — only slightly degraded with respect to corresponding 
values for 0.5 sec delay. At the higher value of gain 

(K£ = 0.0693 in/ft/sec) required to achieve an undamped natural frequency 
of 1 rad/sec, however, the damping ratio is reduced to 0.19 with 1.0 sec 
time delay. Figure 17 shows that the closed-loop damping ratio can be 
restored to about 0.7 with a 1.0 sec delay by reducing the frequency of 


the lead compensation, 1/T^ , to a value corresponding to the heave 
damping 1 /T q . 


All of the foregoing effects and requirements for compensation of time 

delay in height guidance and control systems can be summarized as shown in 

Fig. 18 for a specific value of heave damping 1 /T A = 0.567 rad/sec in the 

°2 

illustrative rotorcraft flight condition. Plotted in Fig. 18 as functions 

of the time delay, t^, are the lead compensation frequency, 1/T L 9 and the 

h 

high frequency gain, -K^Zg required to provide two closed-loop 
characteristic height control frequency bandwidths expressed in terms of 
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Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant A”(s) of Height and Airspeed Control Equations 

While Traveling 

A"(s) = s(s + 1/T Q )(s + 1/T 0 ) - K^[exp(-0.58)](s + 1/T ) • 

1 2 c Hi 
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Figure 1^-. Effects of 0.5 sec Time Delay and Varying Height-to-Collective 
Control Loop Gain with Lead Compensation 1/Tx = 1.13 rad/sec and 
Airspeed-to-Pitch Attitude Control Gain K u = -0.01 rad 0 /ft/sec 
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Real s 


Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant A M (s) of Height and Airspeed Control Equations 

While Traveling 
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Figure 15# Effects of 0.5 sec Time Delay and Varying Height-to-Collective 
Control Loop Gain K£ with Lead Compensation l/T^ - 0.89 rad/sec and 
Airspeed-to-Pitch Attitude Control Gain K u = -0.01 rad 9 c /ft/sec 


50 





Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant A”(s) of Height and Airspeed Control Equations 

While Traveling 
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Figure 1 6. Effects of 1.0 sec Time Delay and Varying Height-to-Collective 
Control Loop Gain with Lead Compensation 1/T. = 1.13 rad/sec and 

Airspeed-to-Pitch Attitude Control Gain K u = -0.01 rad 9 c /ft/sec 
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Figure 17* Effects of 1,0 sec Time Delay and Varying Height-to-Collective 
Control Loop Gain with Lead Compensation l/T^ - 0.56 rad/ sec and 
Airspeed-to-Pitch Attitude Control Gain K - -0.01 rad 9 /ft/sec 






cog, the undamped natural frequency of height regulation, each with closed- 

loop damping ratio, ^ 38 0.7, The two sets of closed-loop characteristics 

are represented by = [0*7;0*5] and [0.7; 1.0]. The effects on 

these closed-loop characteristics of a headwind- to- tail wind shear with 

decreasing altitude will be discussed subsequently. Because all of the 

results in Fig. 18 scale with heave damping, 1/T fl , of the rotorcraft, 

v 2 

they can be normalized by 1/T ^ and presented in the more general form 
depicted in Fig. 19 for the purpose of design. 


Furthermore, since the form of the rotorcraf t's controlled element in 

hover turns is analogous to that for height control, as shown in Table 6, 

the same results, with appropriately revised notation, can be applied to 

hover turn guidance and control with delay as shown in Fig. 20 

(normalized) and Fig. 21 (specialized for effective yaw damping, 

-N' = 0.567 rad/sec). 

r eff 


C. EFFECTS OF HORIZONTAL HIND SHEAR GRADIENT 
(XI HEIGHT CONTROL WHILE TRAVELING IN 
"DOLPHIN" MANEUVERS 


The effects on closed-loop characteristics of a headwind-to-tailwind 
shear gradient, 3u / 3h < 0, with decreasing altitude are shown in Figs. 22 

o 

through 25. The airspeed regulation subsidence, ojJJ = 0.3 rad/sec, 
corresponds with a velocity-to-pitch attitude control high frequency gain 
K y = -0.01 rad 9 c /ft/sec in both sets of figures. Figures 22 and 23 
represent the effects with the lower value of height-to-collective control 
high-frequency gain = 0.0245 in/ft/sec (damping ratio 0.785; undamped 

natural frequency 0.51 rad/sec); Figs. 24 and 25 represent the effects for 
a higher value of = 0.0923 in/ft/sec (damping ratio 0.7; undamped 

natural frequency 1.0 rad/sec). The shear gradient 3u / 3h varies along 

o 

the complex and real root loci in Figs. 22 and 24 as shown in the 

corresponding logarithmic Sigma-Bode root loci, Figs. 23 and 25. The 

results in Figs. 22 through 25 are further expanded in Appendix B to 

represent the effects of two other values of airspeed regulation 

subsidence, a)^ = 0.15 and 0.6 rad/sec, corresponding respectively to two 
other values of velocity-to-pitch attitude control high-frequency gain, 
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NOTATION 


m — = -Z = Heave Damping (rad/sec) 

J-0 2 w 

-Zg c — Collective Control Effectiveness 
ft/(sec 2 -in) or ft/(sec 2 - $) 


Compensation 




5S- r 

C -?. n **£* 


Frequency 


/>» 


/V rr 

* 


Co. ? 




[0.7 • 



fe h^£J ■ [o.7;O.88^t 0 ] 


Height Control Processing Time Delay, T h /T 0 (dmsl) 


Figure 19* Compensatory Height Control High Frequency Gain, 
and Lead Equalization Requirements as Functions of Processing 
Time Delay, Heave Damping, and Collective Control Effectiveness 
Required to Maintain Constant Closed-Loop Characteristics 






NOTATION 



Hover Turn Control 


Processing Delay 
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Figure 20. Compensatory Hover Turn Control High Frequency Gain 
and Lead Equalization Requirements as Functions of Processing 
Time Delay, Yaw Damping, and Yaw Control Effectiveness 
Required to Maintain Constant Closed-Loop Characteristics 
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Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant A”(s) of Height and Airspeed Control Equations 

While Traveling 
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Figure 22. Complex Root Locus of Characteristic Values for the 
Closed-Loop Determinant of Height and Airspeed Control Equations 
While Traveling Showing Effects of Headwind-to-Tailwind Shear 

with Decreasing Altitude 
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Complex Conjugate Roots from Height Regulation 
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Real s 

Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant A”(s) of Height and Airspeed Control Equations 

While Traveling 


A"(s) = [s + 1/T + K u (X a - g)] 

[s 2 + ll/T 9 - Z 5 K^s - (Zj K^/T, )] 
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g v u u or 

3u g /3h < 0, UH-1H at 20 kt, Sea Level, 8000 lb, mid eg, - 0.0923 in/ft/sec 
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-Z u = 0.198 rad/sec 




J c u “0.3 rad/sec 


Figure 24. Complex Root Locus of Characteristic Values for the 
Closed-Loop Determinant of Height and Airspeed Control Equations 
While Traveling Showing Effects of Headwind-to-Tailwind Shear 

with Decreasing Altitude 
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Complex Conjugate Roots from Height Regulation 
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viz., K u = -0.005 and -0*02 rad 9 c /ft/sec. Appendix D provides a 

discussion of airspeed regulation via pitch attitude control and of some 
of the issues affecting the choice of gain K u . 

Not surprisingly, the higher height control bandwidth in Figs. 24 and 
25 provides a higher threshold of instability caused by the wind shear 
gradient. Although the height control bandwidth in Figs. 24 and 25 is 
only doubled with respect to that in Figs. 22 and 23, the wind shear 
gradient required to cause instability is between 2-1/2 and 3 times 
greater in Figs. 24 and 25, depending on the speed regulation subsidence. 

The results from the six examples in Appendix B are summarized in 

Figs. 26 through 28. 

The gradient of the headwind-to-tailwind shear, 3u / 3h <0, is repre- 

s 

sented on the (inverted) ordinate, and the closed-loop speed subsidence 
(rad/sec) in the absence of wind shear is represented on the abscissa 
of Figs. 26 through 28. The closed-loop height control characteristics 

provide damping ratio cjj - 0.785 and undamped natural frequency - 0.51 
in the absence of wind shear in Fig. 26 and Q = 0.7, cd^ = 1.0 without 

wind shear in Fig. 27. Each figure shows a family of curves for the 

relative shear-induced damping ratio and the relative shear-induced 

undamped natural frequency (i.e., height control bandwidth) expressed as a 
percentage of the respective values without wind shear. The shear-induced 
limit of closed-loop stability is represented by the loci for = 0. The 
higher degree of immunity from the compromising effects of the wind shear 
gradient provided by the higher height control bandwidth in Fig. 27 is 
readily apparent. Both Figs. 26 and 27 also demonstrate that a higher 

closed-loop speed subsidence (which provides tighter speed regulation) 

causes the compromised closed-loop height damping ratio to be more 
sensitive to the wind shear gradient than the compromised closed-loop 
height bandwidth; whereas, a lower closed-loop speed subsidence (which 
provides looser speed regulation) causes the compromised closed-loop 
height bandwidth to be more sensitive to the wind shear gradient than the 
compromised closed-loop height damping ratio. Thus an intermediate value 
of the closed-loop speed subsidence will provide a balanced compromise 
among the requirements for speed regulation and immunity from the effects 
of the wind shear on height regulation. 
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for Height Control as a Function of Headwind-to-Tailwind Shear Gradient 
and Closed-Loop Speed Subsidence, if Closed-Loop Height Control Characteristics 
Provide [Sjjaf] = [0.785;0.51 rad/sec) Without Wind Shear 
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and Bandwidth, Whichever is Less, Induced by Wind Shear for Two Different Closed-Loop 

Height Control Characteristics 




To aid in the interpretation of these results and in their use in 
design, we have prepared in Fig. 28 the loci for paired equal percentages 
of relative closed-loop damping ratio and bandwidth, whichever is less, 
induced by the wind shear for the two sets of height control 

characteristics. Two concluding examples will illustrate the use of 
Fig. 28. 

Example 1 : Closed-Loop Height Control Characteristics 

= [0.785;0.51] Without Wind Shear. If one wishes to 

accept no less than 60 percent and 60 percent caused by the 

wind shear, what speed subsidence wJJ should be selected in the 

absence of wind shear and how large a gradient 3u / 3h < 0 can be 

6 

tolerated? The speed subsidence w” = 0.29 rad/sen in the absence 
of wind shear and 3Ug/3h can be at most -0.45/sec 

(-26.6 kt/ 100 ft). 

Example 2 : Closed-Loop Height Control Characteristics 

= [0.7; 1.0] Without Wind Shear. If one wishes to provide 
a degree of height control immunity against a wind shear gradient 
up to 3Ug/3h = -1.0/sec (-59.2 kt/100 ft), what is the largest 
speed subsidence o>^ which will provide balanced relative immunity 
in closed-loop damping ratio cjj and bandwidth o>[|? The speed 

subsidence a)” = 0.55 rad/sec and the balanced relative level of 
height control immunity will be not less than 63 percent ^ and 
63 percent a>^. 

The conclusion at this point is that a closed-loop height-to- 

collective control characteristic having a damping ratio of at least 0.7 
and a bandwidth of at least 0.5 rad/sec will provide insufficient margin 
of immunity against a headwind-to-tailwind shear gradient with decreasing 
altitude in excess of -50 kts/100 ft, although a closed-loop height 
bandwidth of 1.0 rad/sec will provide superior immunity. An intermediate 
range of closed-loop speed subsidence of 0.3 to 0.5 rad/sec in the absence 
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of wind shear provides a more favorably balanced immunity than a high 
subsidence in excess of 0.6 rad/sec, because the wind shear gradient 
3ug/9h < 0 will itself cause the speed subsidence to increase. 
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SECTION IV 


TASK III: AUTOMATIC GUIDANCE FOR AGGRESSIVE NOE MANEUVERS 


A. A TIME-OPTIMAL MODEL FOR AGGRESSIVE 

NOE MANEUVERS 

Four distinctive multiaxis NOE maneuvers, each involving primarily one 
of the four rotorcraft controls as a commanded input, are critical to the 
success of most NOE missions* These maneuvers include the (1) bob-up and 
-down, (2) hover turn, (3) dash-quickstop, and (4) sidestep* Recent Black 
Hawk (UH-60A) flight tests and simulation tests for the purpose of assess- 
ing the fidelity of simulated NOE maneuvers (Refs. 21 through 24) have 
shown that experienced pilots can perform nearly time-optimal bob-ups, 
hover turns, dash-quickstops , and sidesteps when circumstances require 
aggressive maneuvers* Time optimality therefore provides a well-defined 
and validated criterion for corresponding automatically controlled NOE 
maneuvers, if the pilot is to have confidence in the automated maneuvering 
technique* 

Instructions given to the pilots for performing the rapid response 
phase of the bob-ups or -down and hover turns described in Refs. 21 
through 24 were intended to instill a sense of urgency approaching a step 
function for the internal height or heading command while the pilot was 
otherwise engaged in steady-state regulation of his attitude, heading, 
height, and position in the presence of atmospheric turbulence. It has 
been found (Ref. 25 and Appendix A herein) that such combined transient 
and steady-state situations can be represented adequately with a single 
input dual -path structure such as that shown in the simplified diagram of 
Fig, 29, To represent a multiloop situation, the signals shown in this 
block diagram could be considered as vector quantities. The quasi-linear 
steady-state path is the one used for regulating errors caused by random 
inputs or disturbances and illustrated previously in Fig. 10 with the 
symbol Yg. It operates when the error (e) has been reduced within a 
tolerance acceptable to the pilot for the task of regulation. The 
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Rapid 



Figure 29. Structure of the Dual-Mode Model 





feedforward element In Fig. 29 operates on the large transient errors 
Induced by the pilot's internal command or desire to initiate the bob-up 
or -down or hover turn* This feedforward element has been embedded in the 
context of Fig. 10 by the block identified with the symbol Y RR (denoting 
rapid response). 

The roles of the switching and the feedforward element are, in the 
simplest terms, such as to partition the pilot's control strategy into 
three phases, each having a different system organization. As an ele- 
mentary example, consider the typical system step response shown in 
Fig 30. In terms of the three phases, the operation of the dual path 
model can be expressed in the following terms: 

• Transition from quasi-linear path to feedforward path, cor- 
responding to the time delay phase of duration T t « 

• Patterned feedforward response, corresponding to the rapid 
response phase of duration T c . 

• Return to the quasi-linear path, corresponding to the error 
reduction phase of indefinite duration. 

The time delay phase is observable only when the transient forcing func- 
tion is imposed on the pilot unexpectedly from an external source. 

B. HEIGHT CONTROL IN THE BOB-UP OR 
-DOWN RAPID RESPONSE PHASE 


The bob-up or -down rapid response phase begins at the point in time 
that the pilot mentally defines the magnitude of the change in altitude 
required to reach a desired height and decides to displace the collective 
control to begin the maneuver. This is essentially an open-loop command 
which is designed by the pilot to obtain a rate-of-climb or -descent as 
soon as possible while keeping within safe torque limitations and within 
the capability of the rotor speed governor. Initially, no attention is 
given to the problem of stabilizing at the new altitude. 

The most important aspect of the rapid response phase is the pulse- 
like "bang-bang" nature of the control movements, albeit with practical 
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Phase L Rapid Response Phase 

( r t ) ) 


Figure 30 , Typical System Step Response 
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rate limitations* In fact, the pilot^s control displacement (c) is a 
remarkably good approximation to the controller properties of the single- 
input, single-output, time-optimal control system with |c(t)| < M* Here, 
the bound (M) may represent either a physical limit on the control 
deflection or the velocity command or, more likely in the piloted case, an 
implicit restraint imposed by the pilot for the given situation for 
keeping within safe torque and rotorspeed limitations. 

The output of the feedforward element representing the rapid response 
technique for a skilled pilot is peculiar to each controlled element 
form* For the helicopter, the controlled element transfer function, 
Y c (s), representing height response to collective control displacement is 

given with good approximation for our purposes by Y c (s) = K c /s(s + a) 

• 

where K = -Z* t the collective control acceleration effectiveness and a = 
c °c 

-Z w , the heave damping. Ideal time-optimal traces for comparison with the 
actual piloted responses are shown in Figure 31. (Those for Y c (s) ■ K c / s 
and K c /s are from Ref. 25 and Appendix A.) The optimal control pulse 
intervals are unequal in duration for controlled element K c /s(s + a); the 
starting pulse is longer than the final pulse when the damping coefficient 
a is substantial. Suboptimal control techniques by pilots in Refs. 21 
through 24 are evidenced by irregular trapezoidal control, attitude, or 
command pulses of long duration and unequal amplitude. Usually the 
starting pulse is greater in amplitude than the final pulse. Practical 
examples of Y c (s) * K c /s(s + a) and other controlled elements are listed 
in Table 6. Note that the form Y c (s) » K c /s(s + a) applies to the heading 
response-to-pedal inputs in a hover turn with constrained attitude as well 
as to all four maneuvers with the velocity command-displacement hold con- 
trol systems typically required for automatic NOE maneuver control. The 
controlled elements for displacement response to attitude commands at low 
frequencies in the dash-quickstop and the sidestep are given approximately 
by the form Y c (s) = K c /s , because surge damping and sway damping are 
usually very small in rotorcraft at the low speeds which prevail in NOE 
maneuvers. 
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Figure 31% Ideal Time-Optimal Response Characteristics 
(T = Time to Complete Control Program) 



1. Height Control With Vertical Velocity 

Command-Height Hold (VCHH) 

During the rapid response phase of the vertical maneuver with velocity 
command-height hold capability in Fig. 32, the reference for the 
compensatory height regulation loop via the controller is slewed by the 

velocity command integration while the height feedback loop remains closed 
to provide a measure of immunity against disturbances — especially the 
destabilizing effect of a headwind-to-tailwind shear gradient with 
decreasing altitude: 3u g /3h < 0. Concurrently, the three other 

compensatory loops for longitudinal (Fig. 32) and lateral stationkeeping 
(Fig. 33) and heading-hold are closed via the controller throughout the 
rapid response phase of the vertical maneuver. To simplify the subsequent 
analysis, we shall consider only the vertical and longitudinal motions by 
assuming that the controller decouples the lateral-directional motions 

during the vertical maneuver by means of the methods described, for 

example, in Ref. 26. The mathematical details of the following analysis 
appear in Appendix C. Equations C-l through C-4 and the accompanying text 
(pp. C-l through C-3) form prerequisites for interpreting the motion and 
control responses in terras of the symbols in the block diagram of 
Fig. 32. The first subtopic (a) following is expanded in Appendix C, 

pp. C-3 and C-4 using Eqns. C-5 and C-6. 


a. 


Example of design for a special case of VCHH with 
preservation of inherent heave damping, i.e., 1/K^ - 1/T^ 


he 


4 1/T^ - 1/T^ 4 1/Tg^ and no wind shear, i.e. , 3u g /3h = 0 


The analysis in Appendix C shows that the special case of VCHH with 
preservation of heave damping in the absence of wind shear leads to the 
following design constraints 


-Z 



c c 


z s K i "K ' 1/T 

c 


9 


2 
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Figure 32. Block Diagram of Mathematical Models for Covariance Estimation 
in Height Control and Longitudinal Control Tasks 









Measurement Noise 
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Figure 33 . Block Diagram of Mathematical Models for Covariance Estimation 
in Directional Control and Lateral Control Tasks 






















or, in terms of symbols in the block diagram of Fig. 32, 


-h\K "- Z 6 K I h T h" 

c h c c h 


’ * “E * 1/T < 


where 1/T 0 is the inherent heave damping, -Z^ = -Z^ = !/ t q > 

where 1/T 0 ^ t ^ ie inherent heave damping, a function of main rotor pro- 
perties, disc loading, geometry and true airspeed, and Z « , the collective 

°c 

control effectiveness, is also a function of main rotor properties, disc 
loading, geometry, and true airspeed. Equation 05 in Appendix C, p. 03, 
approximates the controlled element form Y c - K c /s(s + a) with K c = a for 
which the generic time-optimal velocity command and rapid height response 
to a step command in height has already been described. Equation 06 in 
Appendix C, p. 04, shows that the required collective control displace- 
ment will be proportional to and in phase with the vertical velocity 
command time history. 


b. Examples of transient responses during a 

time-optimal bob-down for the special Case (a) 
in the absence of wind shear 


For a time-optimal 50 ft bob-down, define amplitude A = 50 ft (in a 
negative sense) and select the dimensionless inverse maneuver urgency 
factor, a 2 A/K c M = 1.0 in Appendix E, Table E-l , p. E-2 , where 
= a = 1 /Tq = 0.567 rad/sec for this example, and M represents the 
absolute magnitude of the vertical velocity command, ± |h c |. Therefore, 
M ~ |^ c | “ A/T 02 = 28.35 ft/sec, and Appendix E, Table E-l, gives the 


descent velocity 


command 


pulse 


time 


interval 


g - 1.61 Tq^ - 2.84 sec, and the final opposite arresting velocity corn- 


starting 
T 

mand pulse time interval T p = 0.511 T 0 - 0.901 sec. Thus the time to 

complete the 50 ft descent will be T = T + T* = 2.12 T Q = 3.74 sec. 

c s r 

Time histories of the ideal optimal velocity command, the change in col- 
lective control displacement from trim, the change in height, and the 
vertical velocity are plotted in Figs. 34 through 37, respectively. The 
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+28.35 ft/sec 


Vertical 

Velocity 

Command 


2.84 see 3.74 sec Tlne (sec) 


Figure 3k. Programmed Time-Optimal Time History of Vertical Velocity 
Command for 50-ft Descent with a Vertical Velocity Command-Height Hold 
Control System Having an Effective Controlled Element Transfer Function 

• 

Y (s) = [h(s)/h (s)] = [0 .567/ s(s + 0.567)] sec 
c c 
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VCHH for the UH-1H; no Shear; 20 kts 


Figure 35. Programmed Optimal Time History of Collective Control 
r °! T o° ft Descent with a Vertical Velocity Command-Height Hold 

Control System Having an Effective Controlled Element Transfer Function 

Y c (s) = [h(s)/h c (s)] = [0.567/s(s + 0.567)] sec 
t 5 coll (s)/h c (s) ] "'TtS’" °- 0610 in/ft/sec 
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Change 

In 

Height 

(«) 



0 12 3^56 

Time (sec) 


7 8 


VCHH for the UH-1H; no Shear; 20 kts 


Figure 36. Time-Optimal Descent of 50 ft with Programmed 
Vertical Velocity Command-Height Hold Control 
for Controlled Element Transfer Function 

• 

Y (s) = [h(s)/h (s)l * [0.567/s(s + 0.567)] sec 
c L c J 
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VCHH for the UH-lH; no Shear; 20 kts 


Figure 37. Phase Plane for Time-Optimal Descent of 50 ft with 
Programmed Vertical Velocity Command-Height Hold Control 
for Controlled Element Transfer Function 

Y c (s) = [h(s)/h c (s)] = [0.567/ s(s + 0.567)] sec 
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2 

required collective displacement authority if -Z. = 9.29 ft/sec -in will 

# 

be ± |S co 2 il = ± | h c | /~Z g Tq = ± 1.73 inches out of a total of 

c 2 

10.7 inches or ± 16.2 percent with respect to the trim displacement which 
is [ ( 13.57—7 .25 ) deg ] /0 . 1 5 deg/percent - 42.1 percent of main rotor 
collective authority for the UH-1H helicopter at 20 kts. 


c. Example of VCHH design Case (a) in a wind shear 
gradient 3u g / 3h = -1.0/sec (-59.2 kt/100 ft) 


Equation C-2 in Appendix C shows that a head- to- tail wind shear 

gradient with decreasing altitude ( 3u g / 3h < 0) in the bob-down will reduce 

the closed-loop characteristic undamped natural frequency a)^, if Z u < 0 , 

which is the usual case for rotorcraft in forward level flight at speeds 

below that for minimum power required*. Figure 38 shows a Sigma-Bode root 

locus of this effect on the pair of characteristic poles in Example a 

above, the special case of VCHH where u)|J * 1 /Tq^ and ^h = * * n the absence 

of wind shear. If 3u / 3h = -1.0/sec, the poles remain stable subsidences, 

but the value of the more critical factor decreases from 

1/Tq - 0.567 rad/sec to 0.12 rad/sec, and the value of the other factor 
°2 

increases from 1 /T a = 0.567 rad/sec to 1.0 rad/sec. 

y 2 

Equation C-4 in Appendix C shows that a head-to-tailwind shear 
gradient with decreasing altitude (3u /3h < 0) in the bob-down will also 

o 

alter the singularities in the numerator of the collective control 
response. Figure 39 shows a Sigma-Bode root locus of the effect on the 
zeroes of this numerator in Example a above. If 3u /3h = -1.0/sec, the 

O 

zero at the origin of the real axis in the complex plane moves to the 
right half plane and remains real at 0.244 rad/sec (factor is -0.244 in 
Fig. 39), and the value of the real zero in the left half plane at 
-1/Tg = -0.567 rad/sec increases (in a negative sense) to -0.811 rad/sec 

(factor is 0.811 in Fig. 39). 


*Reference 27 lists Z u for five helicopters. In forward level flight, 
the largest negative values of Z u occur at 20 kt for four and 30 kt for 
one of the helicopters. 
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Figure 38. Sigma Bode Root I/>cus Showing Effect of Headwind-to-Tailwind Shear with 
Decreasing Altitude: 3Ug/9h > 0 on Characteristic Poles for Vertical Velocity Command-Height Hold 

(VCHH) Control. Longitudinal Position Loop Closed Via Pitch Attitude With Lead Compensation 
1/T x =0.3 rad/sec and High Frequency Gain = -0.02 rad 0 c /ft/sec 
(Complex Conjugate Roots From Stationkeeping Dipole Not Shown at [?" ;oj"] = [0 .7 ; 0 .43 ] 



/l-O- 
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Figure 39* Sigma Bode Root Locus Showing Effect of Headwind-to-Tailwind Shear with 
Decreasing Altitude: 9u^/3h > 0 on Zeros of Numerator for Vertical Velocity Command-Height Hold 

(VCHH) Control* Longitudinal Position Loop Closed Via Pitch Attitude With Lead Compensation 
1/T = 0.3 rad/sec and High Frequency Gain K* = -0.02 rad 0 /ft/sec 

X X c 

(Complex Conjugate Roots From Stationkeeping Dipole Not Shown at = [0.7;0.43] 


d« Examples of transient responses for VCHH 
Case (a) In a wind shear gradient, 

3Ug/ 3h - -1.0/sec (-59.2 kt/100 ft) 

Figures 40, 41, and 42 show the transient responses for collective 
displacement, height displacement, and vertical velocity, respectively, 
for the 50 ft bob-down with VCHH Case (a) In a wind shear gradient 
3Ug/3h = -1.0/sec. The velocity command remains as In Fig. 34 (p. 78). 
The result Is, of course, not time optimal, because the gain of the 
steady-state closed-loop vertical velocity response exceeds unity (it Is 
2.68), and the effective bandwidth of the height-hold loop Is reduced to 
0.12 rad/sec instead of 0.567 rad/sec. Thus the change In height 
overshoots -50 ft during the rapid response phase of 3.74 sec and 
continues to drift downward to approach -130 ft In about 24 sec because of 
the sluggish height-hold loop caused by the wind shear gradient In the 
absence of Integral control during the height error reduction phase. In 
practice, this effect of the wind shear might be counteracted by switching 
to a higher gain compensatory error reduction loop at the end of the rapid 
response phase or by predicting or measuring the wind shear Itself, the 
last two of which alternatives are difficult. Thus It Is preferable to 
redesign the VCHH capability in Case (a) to provide more tolerance to such 
a severe wind shear gradient. 

Example of redesign of VCHH system to 
increase tolerance to a wind shear gradient 
3u g /3h = -1.0/sec (-59.2 kt/100 ft) 

(See Appendix C, pp. C-4 and C-5.) 

f • Redesigned example Case (e) in the absence 
of wind shear ( 3u / 3h = 0) —— 

o 

The analysis for this topic begins on p. C— 6 In Appendix C. In the 
absence of wind shear, the effective controlled element in Eq. C-ll in 
Appendix C, p. C-6 , is still approximately of the form Y c = K,/s(s + a), 



COLL/h c for UH-1H with VCHH; -1.0/sec Headwind-to-Tailwind 

Shear in Descent 


Figure 40. Collective Control Time History with Programmed 
Vertical Velocity Command Appropriate for a Time Optimal Descent 
of 50 ft in the Absence of Wind Shear but Applied in the Presence 
of Wind Shear to VCHH Control Without Integration of Height Error 

in Wind Shear 9u„/8h < 0 


Time (sec) 


Y c ■ h/h c for UH-1H with VCHH; -1 .0/sec Headwind-to-Tailwind 

Shear in Descent 


Figure . Descent Time History of Height with Programmed Vertical 
Velocity Command Appropriate for a Time Optimal Descent of 50 ft 
Showing Stable Drifting Tendency of VCHH Control Without Integration of 
Height Error in Wind Shear 3u /3h < 0 

o 

(Gradient -1.0/sec = -59.2 kt/100 ft) 



Change in Height (ft) 


Y c “ h/ \ for 011-111 w 1 ** 1 VCHH ; ”1 .0/sec Headwind-to-Tailwind 

Shear in Descent 


Figure 42. Phase Plane for Descent with Programmed Vertical 
Velocity Command Appropriate for a Time Optimal Descent of 50 ft 
Showing Stable Drifting Tendency of VCHH Control Without Integration of 
Height Error in Wind Shear 3u /3h < 0 

o 

(Gradient -1.0/sec = -59.2 kt/100 ft) 





where a 4 0.862/sec based on the dimensionless variability* and 
a 4 1.0/sec based on the half-power frequency; but the steady-state gain 
of the velocity response to a velocity command is less than unity (it is 
0.618). 

If the high frequency gain (-Z 5 K£ - 0.567/sec) is increased by a 

c c 

factor 1/0.618 = 1.618 to compensate for the reduced gain in Eq. C-ll, the 
new high frequency gain becomes K c = 1.618 (0.567) = 0.917/sec in 
Eq. C-lla, and if the amplitude and time intervals of the velocity command 
are recomputed to account for the increased K c and increased a in the 
maneuver urgency factor K c M/a 2 A = 1.0 where M= |h c |, then 
|h c | “ a 2 A/K c = (0.862) 2 50/0.917 = 40.51 ft/sec, where a = 0.862 based on 
the dimensionless variability; and T = 1.61/1.0 = 1.61 sec, 
T f = 0.511/1.0 = 0.511 sec, and T £ = 2.12 sec, where a 4 1.0/sec based on 
the half-power frequency. The application of the two different values for 
a is deliberate because of the form of the modified controlled element in 
Eq. C-lla, p. C-7 , Appendix C. 

The modified velocity command (h c ) is shown as a function of time in 
Fig. 43. Transient responses of collective and height displacements for 
this modified (and nearly re-optimized) velocity command are given in 
Fig. 44 and Fig. 45, respectively. Corresponding velocity response is 
presented in the phase plane Fig. 46, which confirms that the height and 
vertical velocity responses have been nearly re-optimized in the absence 
of wind shear. 


*"Bandwidth" is a vague term unless the displayed signal spectrum is 
rectangular. For other spectral shapes, the dimensionless variability can 
be used to define a rectangular bandwidth equivalent, i.e., 


e 


L/ 0 °° * 11 (w)da)] 2 

C du> 


(Ref. 27) 


89 



Time (sec) 


Figure 43. Modified Vertical Velocity Command Representing 
Redesigned Example Case (e) for VCHH in the Absence of Wind Shear 




Time (sec) 


Figure 44. Time History of Change in Collective Displacement 
From Trim in Response to Modified Vertical Velocity Command 
Representing Redesigned Example Case (e) for VCHH in the 
Absence of Wind Shear 


(10.7 inches of collective displacement = 100 percent) 





Time (sec) 


Figure 45. Time History of Change in Height Displacement in Response 
to Modified Vertical Velocity Command Representing 
Redesigned Example Case (e) for VCHH in the Absence of Wind Shear 






g. Conclusion from these examples with VCHH 


The worst case design for the wind shear in Example Case (e) is 
preferable, because it results in a conservative design which can be more 
easily re-optimized as in Example Case (f) by velocity command amplitude 
and time interval adjustments in the absence of wind shear. 

2 . Loss of VCHH Capability 

During the rapid response phase of the vertical maneuver without 
velocity command-height hold capability, the compensatory feedback loop 
via the pilot or controller to the primary rotorcraft control, in this 
instance, the collective control, is open. It will therefore be of 
comparative interest subsequently to examine the effect of a headwind-to- 
tailwind shear with decreasing height on the time-optimal bob-down during 
the rapid response phase with open height control loop. This situation 
would apply in the event of a passive open failure of the height feedback 
loop in the vertical velocity command-position hold control system, 
Fig. 32 (p. 75). The three other compensatory loops for longitudinal 

(Fig. 32) and lateral stationkeeping (Fig. 33, p. 76) and heading-hold 
are, however, presumed to be closed via the controller (human or 
automatic) throughout the rapid response phase of the bob-down maneuver. 
To simply the subsequent analysis, we shall restrict it to consider only 
the vertical and longitudinal motions by assuming that the controller 
decouples the lateral-directional motions perfectly. 

a. Longitudinal Hovering Position Control in Vertical 

Unmask and Remask Maneuvers with Open-Loop Height Control 

Figure 47 shows the effects while maintaining a longitudinal hovering 

position, i.e., stationkeeping, with open-loop height control, of varying 

longitudinal posi t ion-to-pi t ch attitude control high-frequency gain, K£, 

with constant lead compensation, 1/Ty =0.3 rad/sec on the closed-loop 

x 

dynamic characteristics of longitudinal translation. Three values of the 
gain and their corresponding closed-loop damping ratio and undamped 
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Figure 47 . Effects of Varying Longitudinal Position-to-Pitch Attitude 

Control Loop Gain with Lead Compensation 1/T, *0.3 rad/sec 

and Without Height Regulation 


natural frequencies are identified. There is no time delay in the mathe- 
matical model represented in Fig. 47. Likewise, the higher frequency 
pitch attitude dynamics vested in the residual oscillatory "pendulum com- 
ponent" characterized by the u> u -u)' dipole (refer to Table 6, p. 38, and 

^ * 

Appendix D, p. D-4) are not represented in Fig. 47 . The predominant char- 
acteristic dynamic mode of longitudinal translation will be of second 
order, because the heave damping subsidence is suppressed by the transla- 
tion zero, 1/T u = 1/T 0 . 

The higher frequency attitude dynamics associated with the o> u - 

dipole shown in Table 6 and discussed in Appendix D are nevertheless im- 
portant, because the amplitude of the residual pendulum oscillation tends 
to increase and eventually to become objectionable to the pilot as the 
gain increases beyond the range shown in Fig. 47. Criteria for op- 
timizing the gain K^, and its accompanying dynamic effects on translation 

and attitude are discussed in Refs. 29 through 32 in the context of pilot 
opinion ratings derived from simulation and flight test experiments. 

b. Transient Responses During a Time Optimal 
Bob-Down Without Wind Shear While 
Stationkeeping with Open-Loop Height Control 

A more complete set of time histories and phase plane portraits during 
the rapid response phase of a time optimal 50-ft bob-down while station- 
keeping in a steady 20 kt head wind (without any shear) is shown in 
Figs. 48 through 54. There is no VCHH capability. The heave damping 


*"Higher frequency” refers to the frequency in the complex zero of the 
"pendulum dipole" represented by 


0) 



where g is gravitational acceleration and, if S B is expressed in terms of 
swashplate angle in radians, X 
of NOE operations. 


g = g in hovering and at low speeds typical 
B 
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2.84- sec 3.74 sec 


Time (sec) 


- h/COLL for the UH-1H; no Shear; 20 kts 


Figure 48. Programmed Optimal Time History of Collective Control 
for 50 — ft Descent Controlled Element Transfer Function 

Y c (s) = [h(s)/6 (s)] - [9.29/s(s + 0.567)] ft /in 





Time (sec) 


Y c - h/COLL for the UH-1H; no Shear; 20 kts 


Figure ^9. Time-Optimal Descent of 50 ft with Programmed Open-Loop 
Collective Control for Controlled Element Transfer Function 

Y c (s) = [ h ( s ) /6 C0LIj ( s ) ] = [9.29/s(s + 0.567)] ft/in 
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Y c = h/COLL for the UH— 1H; no Shear; 20 kts 


Figure 50 . Phase Plane for Time-Optimal Descent of 50 ft with 
Programmed Open-Loop Collective Control for Controlled 
Element Transfer Function 

Y c (s) = [h(s)/6 C0LL (s)] = [9.29/s(s + 0.567)] ft/in 
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UH-1H at 20 kts, Sea Level, 8000 lbs, mid c.g., no shear. 
Closed-Loop Transfer Function of Controlled Element for 
* -0.02 rad 9 c /ft/sec and 1/T^ «■ 0.3 rad/sec Without Wind Shear: 

N X 

<$ C0LL 0.01126(s + 0.3)(s - 0.0600) rad 

A” (s + 0.543) [ s 2 + 2(0.715)(_0.449)s + (0.449) 2 ] in 


Figure 51» Change in Commanded Pitch Attitude From Trim During 
Stationkeeping Accompanying Time Optimal Descent of 50 ft with 
Programmed Open-Loop Collective Control 
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3 deg/sec 
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Figure 52. Phase Plane for Change in Commanded Pitch Attitude From Trim During Stationkeeping 
Accompanying Time-Optimal Descent of 50 ft with Programmed Open-Loop Collective Control. 

UH-1H at 20 kts, Sea Level, 8000 lbs. No Shear. 
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UH-1H at 20 kts, Sea Level, 8000 lbs, mid c.g., no shear. 
Closed-Loop Transfer Function of Controlled Element for 
» -0.02 rad 9 c /ft/sec and 1/T^ =0.3 rad/sec Without Wind Shear: 

N *C0H 0 .563s(s - 0.0600) ^ 

A ” s(s + 0.543) [ s 2 + 2[0.715)(0.449)s + (0.449) 2 ] in 

Figure 53* Change in Regulated Longitudinal Displacement From 
Stationkeeping Position Accompanying Time Optimal Descent of 50 ft with 
Programmed Open-Loop Collective Control 
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Time, 1 6 sec 


Aft 



(-Z w ) in the example is 0.567 rad/sec, and the collective control accel- 
eration effectiveness (-Zg ) is 9.29 ft/sec -in. The maneuver is 

performed by the open-loop time optimal collective control pulse sequence 
in Fig. 48, using ± 1.73 inches of collective authority from trim 
(10.7 inches * 100%). Changes in height and vertical velocity are dis- 
played in Figs. 49 and 50, pitch attitude and rate command required for 
stationkeeping in Figs. 51 and 52, and transient longitudinal displacement 
deviation and velocity during stationkeeping in Figs. 53 and 54. 

Figures 49 and 50 show that the 50 ft descent is completed in 
3.74 sec, with a residual velocity of about -2 ft/sec, which decays to 
zero at 6 sec even with the compensatory height regulation loop open. The 
final unregulated height is -54 ft. The peak descent velocity at the 
collective control reversal time (2.84 sec) is -33.5 ft/sec. 

Figures 51 and 52 show that the change in commanded pitch attitude 
required for automatic longitudinal stationkeeping regulation is within 
±0.9 deg and returns to trimmed attitude within 10 sec. The extreme 
pitch rates are -1.1 deg/sec (nose down) initially and 2.5 deg/sec (nose 
up) at the time of collective control reversal (2.84 sec). 

Figures 53 and 54 show that the extreme values of displacement from 
the regulated longitudinal stationkeeping position are -1.1 ft (aft) at 
3 sec and 0,7 ft (forward) at 7 sec. The initial stationkeeping position 
is restored in 16 sec, and the extreme values of regulated longitudinal 
velocity are -0.5 ft/sec (aft) at 1.5 sec and 0,9 ft/sec (forward) at 
3.74 sec, when the descent is completed. 

c. Effect of Headwind-to-Tailwind Shear with Decreasing 

Altitude on Longitudinal Stationkeeping During Bob-Down 

While Stationkeeping with Open-Loop Height Control 

(1) Complex and Bode Root Loci of Closed-Loop Characteristics. 
Starting with each of the closed-loop characteristics represented 
in Fig. 47 and Figs. 55 and 56 then show the effect of a head- 
wind-to-tailwind shear with decreasing altitude on the closed- 
loop stability of the longitudinal position control loop during 
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Figure 55. Complex Root Locus Showing Effect of Headwind-to-Tal lwind Shear 
with Decreasing Altitude on Longitudinal Stationkeeping During Bob-Down 
With Open-Loop Height Control: 3Ug/3h < 0. Longitudinal Position Loop Closed with 

Lead Compensation 1/T L = 0.3 rad/sec. UH-1H at 20 kts, Sea Level, 8000 lb, mid eg. 



Complex Conjugate Roots from Stationkeeping Dipole. 
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Figure 56. Sigma-Bode Root Locus Showing Effect of Headwind-to-Tailwind Shear with Decreasing 
Altitude on Longitudinal Stationkeeping During Bob-Down with Open-Loop Height Control: 

3Ug/ 3h < 0. Longitudinal Position Loop Closed Via Pitch Attitude with Lead Compensation 
1/T^ =0.3 rad/sec and High Frequency Gain = -0.02 rad 0 c /ft/sec. 

UH-1H at 20 kts, Sea Level, 8000 lb, mid eg. 



the rapid response phase with open-loop height control. The 
shear gradient causes immediate instability when there is no 
height regulation loop during the bob-down. Values of the diver- 
gence poles are shown along the real axis as a function of the 
gust gradient 3Ug/ 9h for a particular value of high-frequency 
gain = -0.02 rad/ft/sec. For example, if the shear gradient 
3ug/ 3h is -0.4/sec (-23.7 kts/100 ft), the divergence pole will 
be 0.12 rad/sec, which has a time to double amplitude of about 
5.8 sec. If this shear gradient increases to -1.0/sec 

(-59.2 kt/100 ft), the divergence pole will roughly double to 
0.24 rad/sec, which has a time to double amplitude of 2.9 sec. 
The results in Figs. 55 and 56 are not sensitive to the particu- 
lar value of high-frequency position loop gain K^. Likewise, the 
results in Figs. 55 and 56 are not sensitive to the effects of 
delay in the longitudinal position loop. In fact, the results in 

Figs. 55 and 56 are substantially the same as would be obtained 

with only the attitude loops closed (Ref. 33). 

(2) Transient Responses During a Time-Optimal Bob-Down in a Headwind- 
t o-Tailwind Shear with Decreasing Altitude. Time histories and 
phase plane portraits of the divergent bob-down and station- 

keeping maneuvers with a shear gradient 3u /3h = -1.0/sec are 

o 

shown in Figs. 57 through 63 for comparison with Figs. 48 through 
54. Since the 50-ft bob-down requires only 3.74 sec, the effect 
of the relatively rapid divergence can be nulled by closure of 
the compensatory height loop at the end of the maneuver. None- 
theless, the previous conclusion from the examples with velocity 
command-height hold is reinforced, viz., that the designer should 
try to provide reliable velocity command-position hold capability 
in all axes for superior immunity to wind shear in an automatic 
guidance and control system. 
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Y c = h/coll for the UtT-lH; -1.0/sec headwind-to-tailwind shear 
in descent (-59.2 kt/100 ft shear gradient) 

Controlled Element Transfer Function 

Y ( s) = [h(s)/5 coU (s)] = [ 9 . 29 / ( s - 0.243)(s + 0.793)] ft/in 


Figure 57, Programmed Optimal Time History of Collective Control 
Appropriate for 50-ft Descent in the Absence of Wind Shear 
but Applied in the Presence of Wind Shear 
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Y c * h/coll for the UH-1H; -1.0/sec headwind-to-tailwind shear 
in descent (—59 • 2 kt/100 ft shear gradient) 


Controlled Element Transfer Function 
Y c (s) = [h(s)/6 con (s)] = [9.29/(s - 0.243)(s + 0.793)] ft/in 


Figure 58. Descent Time History with Programmed Open-Loop Collective 
Control Appropriate for a Time Optimal Descent of 50 ft Showing 
Destabilizing Effect of Wind Shear 3u /9h < 0 
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Y = h/coll for the UH-1H; -1.0/sec headwind-to-t ailwind shear 
c 

in descent 

Controlled Element Transfer Function 

Y (s) = [h(s)/5 coll (s) ] = [9.29/(s - 0.243)(s + 0.793)] ft/in 

Figure 59. Phase Plane for Descent with Programmed Open-Loop Collective 
Control Appropriate for a Time Optimal Descent of 50 ft Showing 
Destabilizing Effect of Wind Shear 3Ug/3h < 0 

(Gradient -1.0/sec = -59.2 kt/100 ft) 
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UH-1H at 20 kts, Sea Level, 8000 lbs, -1.0/sec headwind-to-tailwind shear 
in descent (-59.2 kt/100 ft shear gradient) 

Closed-Loop Transfer Function of Controlled Element for 
K- = -0.02 rad 9/ft /sec and 1/T^ =0.3 rad/sec Without Wind Shear 


6 coll ... [0.01126(8 + 0.3)(s - 0.4624)(s + 0.5224)] 

a" “ 5 — r rad/in 

(s + 0.243)(s + 0.7934)[s + 2(0.7075)(0.4336)s + (0.4336) Z 


Figure 60. Change in Commanded Pitch Attitude From Trim During 
Stationkeeping with Programmed Open-Loop Collective Control Appropriate 
for a Time Optimal Descent of 50 ft Showing Destabilizing Effect of 

Wind Shear 3Ug/3h < 0 


111 


Nose Up 


e 

<y 


° g> 04 
ft ON 


M I O 


m 

o 

a 

in 

rP 

w 

4-» 

•H 

CM 

m 

o 

u 

u 

CO 

OJ 

pC 

T3 

0) 

O 

o 


I 

T—f 

o 

CD 

CO 

+ 


ft 

CO 

VP 

ft 


CO 

CM 

'S 

T3 

4-1 

H 

CO 

P 


O 

o 

-p 

•H 

-P 

c 

•H 

N 

O 

o 

ro 

<r 

CN 

vD 


-P 

< 

i— ( 

•H 

CO 

MP 

O 

o 

n 

<r 

• 

o 


O 

-P 

4-> 

c 

i 

o 

ft 

1 

o 

4-> 

o 

•H 

4—1 

X 

hJ 

E-h 

i 

CO 

o 

<D 

1 

o 

c 


v ~' 


TJ 

C! 


m 


<y 1 O 

fl &0 0> 

C bO O C m 

3 ,£ C .H *H \ 

S O <D *rJ -P ft TJ 
g +> -p p cd 0) co 
O -H d pJ -P 0) u 

o ft ft Q w ^ 


112 


Figure 61. Phase Plane for Change in Commanded Pitch Attitude From Trim During, Stationkeeping with 
Programmed Open-Loop Collective Control Appropriate for a Time Optimal Descent of 50 ft 
Showing Destabilizing Effect of Wind Shear 9u 0 /3h < 0 


Time, 3-74 see Time, 6 sec 
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Figure 62. Phase Plane for Regulated Longitudinal Displacement from Stationkeeping Position with 
Programmed Open-Loop Collective Control Appropriate for a Time Optimal Descent of 50 ft 
Showing Destabilizing Effect of Wind Shear 3u /3h < 0 






UH-1H at 20 kts, Sea Level, 8000 lbs, -1.0/sec headwind-to-tailwind shear 
in descent (-59.2 kt/100 ft shear gradient) 

Closed-Loop Transfer Function of Controlled Element for 
K* = -0.02 rad 9/ft/sec and 1/T L = 0.3 rad/sec Without Wind Shear 
x x 


N 


coll 


0.563(s ~ 0.4624)(s - 0.5224) 

(s + 0.243)(s + 0.7934) [ S 2 + 2(0.7075)(0.4336)s + (0.4336) 2 


rad/ in 


Figure 63. Change in Regulated Longitudinal Displacement from 
Stationkeeping Position with Programmed Open-Loop Collective Control 
Appropriate for a Time Optimal Descent of 50 ft Showing Destabilizing 

Effect of Wind Shear 9Ug/3h < 0 
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C. EFFECT OF HEADWIND-TO-TAILWIND SHEAR WITH FORWARD 

DISPLACEMENT ON LONGITUDINAL POSITION CONTROL 

Figures 64 and 65 show the effect of the headwind-to-tailwind shear 
with forward displacement, 3u g /9x > 0, on the closed-loop dynamics of the 
longitudinal position— to— pitch attitude control loop* The lead compensa- 
tion 1 / t l x is constant at 0.3 rad/sec, and the gain 
K x = 0*01. rad 9 c /ft/sec. Note that this is a horizontal displacement- 
dependent shear; consequently, the longitudinal position control loop 
tends to resist the gradient of gust velocity with position. Since the 
horizontal-displaceraent-dependent shear acts on the vehicle through surge 
damping (X u ), the very low value of surge damping for this particular 
example makes the closed— loop characteristics particularly resistant to 
the effects of this type of shear. Although this shear causes an im— 

®®^^®te divergence, the divergence pole is so small , even for extremely 
^■®^8® shear gradients, that the time to double amplitude is very long. 
Even if the surge damping were doubled, half of the gradient shown in 
Figs. 64 and 65 would produce the same effects on the closed-loop sys- 
tem. Even then, the smallest of the shears is so extraordinarily large as 
to be extremely rarely encountered. 

The conclusion, then, for longitudinal position— to— pitch attitude 
control is that it provides a virtual immunity to horizontal-displacement- 
dependent wind shear effects. 

This completes our summary of activities on Task III. We shall com- 
plete this discussion of all Phase I activities with a summary of findings 
in Section V for each of the tasks that constitute this phase of the 
study. 
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P30.5/S 


Su g /Sx 


15.7/sec 

I 30 . 5 /s 


UH-1H at 20 kt, Sea Level, 8000 lb. Mid eg 


Longitudinal Posi tion-to-Pitch Attitude Control Closed with Lead 
Compensation 1/T^ = 0.3 rad/sec and Gain - -0.01 rad 9 c /ft/sec 


Figure 64. Complex Root Locus Showing Effect of Headwind-to-Tailwind 
Shear with Forward Displacement on Longitudinal Position 
Control: 9u / 8x > 0 



Complex Conjugate Roots from Stationkeeping 
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Figure 65. Sigma-Bode Root Loci Showing Effect of Headwind-to-Tailwind Shea 
with Forward Displacement on Longitudinal Position Control: 3u„/3x > 0 



SECTION V 


SUMMARY OF FINDINGS 


This section summarizes the principal findings for each of the tasks 
that constitute this phase of the study. 

A. TASK Is STORAGE AND UPDATING OF TERRAIN PROFILES 

The study investigated the representation of the guidance command 
vector for a planned flight profile in terms of a combination of poly- 
nomial and harmonic functions of along-course and across-course 
coordinates, i.e., a vector of truncated Taylor and Fourier series. Ex- 
amples of flight profiles were sampled from a 5 km x 5 km digitized 
terrain model having 1 dm resolution. The examples for which finite 
Fourier transforms were obtained exhibited effective spatial bandwidths of 
about 0.2 rad/dm. At 20 kt forward speed, spatial bandwidth in rad/dm and 
temporal bandwidth in rad/sec are virtually identical. Thus we can expect 
a typical guidance and control system to follow an offset profile over 
these samples with good precision of 20 to 40 kt forward speed. 

The techniques examined for representing a planned flight profile 
using Fourier series descriptors achieved data compression ratios between 
8:1 and 10:1 for the samples of vertical profiles from the particular 
terrain data base used. The resulting recommendations for compressing 
flight profile data storage requirements are summarized in Table 5, p. 27. 

A procedure is outlined for updating and modifying the stored guidance 
command flight profiles in real time to provide unexpected sensed obstacle 
and threat avoidance using constrained time-optimal maneuver strategy. It 
remains to develop in the second phase of this study a rationale and pro- 
cedure for correcting potential conflicts caused by unexpected obstacles 
and threats and returning to the flight plan safely by means of transition 
logic between the anticipative trajectory coupler and the constrained 
time-optimal maneuvers for avoiding obstacles. 
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B. TASK II: AUTOMATIC GUIDANCE FOR FOLLOWING FLIGHT PROFILES 

The automatic guidance algorithms will employ estimates of anticipated 
» lateral, and longitudinal deviation commands required to follow 
the desired path, course, and schedule, respectively, with feedforward 
compensation. Preview of both stored and sensed terrain and obstacle data 
is therefore essential to compensate for the inevitable and predictable 
processing delays and rotorcraft response lags. 

That part of the guidance error associated with following commands may 
be reduced by an automatic pursuit feedforward guidance adjustment anal- 
ogous to that which the pilot will adopt in following a fluctuating curved 
course if sufficient visual preview of the source slope and curvature are 
available. If incorporated in automatic guidance, this practical pursuit 
adjustment will enhance pilot acceptance of automatic guidance in follow- 
ing nap-of-the-earth profiles with precision. A stored data base for the 
flight profile will provide a practical basis for the preview necessary to 
implement this pursuit adjustment* 

A significant degree of immunity from wind shears and turbulence can 
be provided by automatic velocity command-position hold guidance and con- 
trol systems. A closed-loop height response-to-collective control 
bandwidth of 1.0 rad/sec with 0.7 damping ratio will provide superior 
immunity against a headwind-to-tailwind shear gradient with decreasing 
altitude in excess of -50 kt/100 ft, although a bandwidth of only 

0.5 rad/sec will be fairly effective in countering lesser gradients. An 
intermediate range of closed-loop speed subsidence of 0.3 to 0.5 rad/sec 
will help to provide a more favorably balanced immunity than a higher 
subsidence in excess of 0.6 rad/sec, because the wind shear gradient 
3iig/ 3h < 0 will itself cause the speed subsidence to increase. 

C. TASK III: AUTOMATIC GUIDANCE FOR 

AGGRESSIVE NAP-OF-THE-EARTH MANEUVERS 

Four distinctive multiaxis NOE maneuvers, each involving primarily one 
of the four rotorcraft controls as a commanded input, are critical to the 
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success of most NOE missions. These maneuvers Include the (1) bob-up and 
-down, (2) hover turn, (3) dash-quickstop, and (4) sidestep. Recent Black 
Hawk (UH-60A) flight tests and simulation tests for the purpose of assess- 
ing the fidelity of simulated NOE maneuvers have shown that experienced 
pilots can perform nearly time-optimal bob-ups, hover turns, dash- 
quickstops, and sidesteps when circumstances require aggressive 
maneuvers. Time optimality therefore provides a well-defined and vali- 
dated criterion for corresponding automatically controlled NOE maneuvers, 
if the pilot is to have confidence in the automated maneuvering technique. 

The output of the feedforward element representing the rapid response 
technique for a skilled pilot is peculiar to each controlled element 
transfer function form Y c (s). The form Y c (s) = K c /s(s+a) applies to the 
heading response-to-yaw control inputs in a hover turn with constrained 
attitude as well as to all four maneuvers with the velocity command- 
displacement hold guidance and control systems typically required for 
automatic NOE maneuvers. The optimal control pulse intervals are unequal 
in duration for performing aggressive maneuvers with controlled elements 
having the form K c /s(s+a); the starting pulse is longer than the final 
pulse when the damping coefficient a is substantial. Necessary con- 
straints that contribute to suboptimality are rotor torque limitations and 
rotor speed governor recovery capability, control displacement and rate 
limitations, attitude displacment and rate limitations, and acceleration 
limitations in uncoordinated maneuvers. 

Among the examples of design for time-optimal velocity command-height 
hold in a wind shear environment, the "worst case" design to preserve 
heave damping in an extreme wind shear 3Ug/3h - -59 kt/100 ft is prefer- 
able, because it results in a conservative design that can be more easily 
re-optimized by velocity command amplitude and time interval adjustments 
in lesser gradients or in the absence of wind shear. Thus the designer 
should try to provide reliable velocity command-position hold capability 
in all axes for superior immunity to wind shear in an automatic guidance 
and control system for NOE operations. 

The longitudinal position control loop tends to resist the gradient of 
longitudinal gust velocity with horizontal displacement, 3u /3x. Since 

O 
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the horizontal -displacement-dependent shear acts on the vehicle through 
surge damping (X u ), the very low value of surge damping for rotorcraft 
makes the closed-loop characteristics particularly resistant to the ef- 
fects of this type of shear. Although this type of shear causes an 
immediate divergence, the divergence pole is so small, even for extremely 
large shear gradients, that the time to double amplitude is very long. 

The conclusion, then, for longitudinal position-to-pitch attitude 
control is that it provides a virtual immunity to horizontal-displacement- 
dependent wind shear effects. 


This completes our summary of findings on Tasks I, II, and III. We 
shall conclude this discussion of all Phase I activities with a summary in 
Section VI of pilot— centered considerations for monitoring automatically 
controlled NOE flight. 
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SECTION VI 


CONCLUDING SUMMARY OF PILOT-CENTERED CONSIDERATIONS FOR MONITORING 
AUTOMATICALLY CONTROLLED NOE FLIGHT 


A. AUTOMATIC-MANUAL CONTROL RESPONSE COMPATIBILITY 

Automatic pursuit guidance for NOE operations will improve pilot 
acceptance of automatic course-and-prof ile following, because pursuit 
guidance is compatible with manual piloting technique if sufficient 
preview of the course-and-prof ile is available* Automatic guidance of 
aggressive tactical maneuvers for NOE operations based on principles of 
time-optimality (with appropriate control rate, attitude rate, and atti- 
tude constraints) will improve pilot acceptance, because time-optimality 
is compatible with measured manual piloting techniques under visual condi- 
tions in flight tests. Automatic control algorithms should weight the 
various motion and position feedback and feedforward signals so that the 
resulting rotorcraft motions are not disharmonious or dangerous (e.g., no 
excessive attitude overshoot to correct velocity or position error) but 
are similar to those experienced under visual manually controlled 
conditions. It also implies that a flight director be compatible with the 

automatic system. The control techniques for most NOE operations will be 

. 

appropriate for speeds below that for minimum power required, viz., h, h 
5 C ; x,x, 9 -*• 5 b ; y, y, 4 ■*> S A ; i|>, > 5 p 

B. STATUS DISPLAYS 

Status displays should provide cues similar to those in visual flight 
and should be compatible with any command display to inspire pilot confi- 
dence in an automatic control system. In particular for NOE operations, 
pursuit displays such as the pole-track (Refs. 34 through 36), channel 
(Refs. 37 and 38), and their antecedents (e.g., see the review in Ref. 39) 
will become necessary for course-and-prof ile following and maneuvering. 
Pursuit displays will very likely have to be complemented with auditory 
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command and status information during automatically guided constrained 
time-optimal maneuvering. 

C. DECOUPLED CONTROL RESPONSES 

Rotorcraft responses to control inputs are inherently coupled when 
there is but a single main rotor. The predominant coupled responses are 
those in yaw, pitch, and roll due to collective control inputs. Depending 
on the airspeed and the ambient winds, the coupled heading and attitude 
changes can also be accompanied by unwanted translational motions. The 
automatic system should compensate for collective coupling by 

simultaneously applying pedal, lateral stick, and longitudinal stick 
inputs when moving the collective control. 

The prevalent direct and cross-coupled responses to rotorcraft con- 
trols are listed in Table 7 for a single main rotor. The columns are as 
follows: Column 1 identifies the name of the control, Column 2 identifies 

the direct response. Column 3 identifies the predominant cross— coupled 
responses, and Column 4 identifies the crossfeed required to decouple the 
response. Compensation for the cross— coupled responses to control inputs 
may consist of either or both crossfeeds between the controllers and/or 
feedbacks of the responses to the controllers. Feedback offers the 
additional advantage of compensation for responses to unmeasurable distur- 
bances. Feedback also alters the automatic— manual control response 
compatibility of the vehicle. Using crossfeeds between the controllers 
reduces the unwanted vehicle response before it has time to build up. 
Thus the hybrid approach to using both crossfeeds and feedbacks offers a 
cost effective compromise. 

This concludes our recitation of Phase I results. We have demon- 
strated initially the technical feasibility of the objectives of this 
research in terms of the specific findings of each task in Phase I. 
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TABLE 7. COMPENSATION FOR CROSS -COUPLED RESPONSE TO CONTROLS 
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APPENDIX A 


MODELS FOR TRANSIENT INPUTS 
(Excerpts from Ref. A-l) 


A. INTRODUCTION 


The simplest nonstationary control situation is one in which a highly 
trained, but nonalerted, subject operating a constant-coefficient linear 
controlled element in a compensatory system is confronted with a randomly 
occurring step input. In spite of the simple circumstances, the overall 
behavior is complicated when a variety of controlled elements is con- 
sidered. The skilled operator's output is peculiar to each controlled 
element form. The system response is, however, less variant in that it 
tends to duplicate, after a time delay, the forcing function. Thus, the 
system output to a unit step forcing function shown in Fig. A-l is typi- 
cal. This operator response can be analyzed by considering the three 

phases separately. When described in terms of a block diagram that re- 
lates stimulus to response, each temporal phase can be conceived of as 
having a different system organization. The block diagram structure indi- 
cates the dynamics of the association between the pilot's response and the 

actual or effective inputs that he is operating on to generate that 

response • 

For tracking random inputs with occasional step inputs, a dual-mode 
model of the operator is appropriate. The basic structure of the model is 
given in Fig. A-2. The quasil inear steady-state path is the one used for 
tracking random inputs when the error e is less than or equal to an error 
threshold e t , while the feedforward element operates on the random- 

occurring step inputs when the error e is greater than the threshold e t . 
The basic structure thus incorporates mode switches for the two path- 
ways. In terms of the three temporal phases, the successive action 

structures of the dual— mode model are! 
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Figure A-l. Typical System Step Response 
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Figure A-2* Structure of the Dual-Mode Model 





• Transition from quasilinear mode to feedforward mode, corres- 
ponding to the time delay phase 

• Patterned feedforward response, corresponding to the rapid 
response phase 

• Quasilinear mode, corresponding to the error reduction phase. 

The duration of the time delay phase has some minimum value, and its 
unimodal distribution is therefore skewed toward the larger values. The 
time delay, T t , is generally longer than the steady-state effective time 
delay, t e , in tracking. In particular, 

x t = T e + x d ( A— 1 ) 

where T t is the time delay phase duration, x e is the effective time delay 
in steady-state tracking, and r d is the decision time. During the deci- 
sion time, the operator makes the pertinent decisions regarding the shape 
and magnitudes of the feedforward response. At the end of the time delay 
phase, the feedforward element generates the proper response to the step 
input command, giving the rapid response phase. At the end of this per- 
iod, the error is small, and the operator switches to a quasilinear 
tracking mode in the error reduction phase. 

Much is known about the quasilinear controller; therefore, attention 
here will be confined to the parallel feedforward path employing the rapid 
response algorithm used in response to step inputs. 

B. A DUAL-MODE CONTROLLER MODEL 

An important aspect, obvious from the step response data in Refs. A-l 
through A-5, is the bang-bang nature of the stick deflection control move- 
ments. This property leads us to the pertinent problem of optimality of 
the operator and his related performance indices. One explanation is that 
the operator is optimal or suboptimal relative to the minimum time criter- 
ion. Consider an nth order single input single output control system with 
| c ( t ) | < M, where the scalar M may represent either a physical limit on 
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the stick deflection or, more likely, an implicit restraint imposed by the 
operator for the given situation. In any case, it represents a magnitude 
constraint on the control input. For i(t) ■ constant, the time optimal 
control has the following properties: 


• The control c(t) is bang-bang, i.e., c(t) = +M or -M. 

• There are at most (n-1) switchings (i.e., +M or -M or vice 
versa) for systems with n real eigenvalues. 

• The switching logic is dependent on the order of the con- 
trolled element. In general, the switching surface is a 
nonlinear function of the state variables. 

• For a given initial condition of the state variables, there is 
one unique control c(t). 


For the problem at hand, there is a specific type of initial condition 
of the system state vector, namely: 

e(0) = input height 
e(0) = 0 
e(0) = 0 


The terminal state is the origin. 

In order to characterize the degree to which the available step re- 
sponse data with c(t) = ±M is time optimal, certain invariance conditions, 
one for each controlled element, are obtained by solving a two-point boun- 
dary value problem. These are stated and described in Table A-l without 
presenting their derivation; thus, let 

T c =* time to complete the force response (i.e., duration of 
the stick response correction for step inputs) 

M * average absolute amplitude (for each Y ) of the stick 

response assuming it to be bang-bang witn equal positive 
and negative amplitudes 
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TABLE A-l 


INVARIANCE CONDITIONS FOR TIME OPTIMALITY WITH A STEP INPUT 



.-6 






A = amplitude of step input 

K c = controlled element gain 

(*) = time optimal value of the parameter in parenthesis 

o 

Assuming a wide band neuromuscular system (or controller) response, 
the ideal time-optimal step response character for differing controlled 
elements is shown in Fig. A-3. Note that the smoothing effect of the 
neuromuscular system (or controller) would round off the corners in the 
Fig. A-3 responses. Note further that the control movement starts after 
the end of the delay time phase, and time optimality pertains to that 
period of control only. Refs. A-l through A-5 present a comparison of the 
actual data to the optimal. On the basis of the comparisons in the cited 
references, it may be safely concluded that the step response behavior of 
operators is nearly time optimal. 

C. SUMMARY 


In conclusion, a time-optimal control model is one possible idealiza- 
tion for the feedforward step response path of the dual-switched-raode 
model for the operator. The complete model is presented in Fig. A-4 , as 
one explanation of available data. The quasilinear path is the usual 
operator describing function for compensatory steady-state tracking of 
random inputs. The feedforward parallel path represents the control plus 
decision model of the operator in response to step inputs. The nonlinear 
error sensing blocks in Fig. A-4 automatically route the error signal 
through the appropriate channel based upon whether e \ e T (e T is some 
threshold magnitude of error unique to the task). 

The control logic for each different controlled element, and as a 
function of the error state e_ [e = col (e, e, ...)], is given in Table A-2 
for time-optimal response. Note that M, the constraint on the control 
input, is some function of the step input height, controlled element gain, 
and its order. 
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Figure A-3. Ideal Time -Optimal Response Characteristics 
(T = Time to Complete Control Program) 
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Figure A-4 .DUAL-MODE CONTROLLER MODEL 






TABLE A- 2 


CONTROL LOGIC FOR VARIOUS CONTROLLED ELEMENTS 


CONTROLLED ELEMENT 
Y c 

CONTROL LOGIC 
f(e) 

K c 

(A/MK c )e(t) 

K c /s 

e(t) 

K c /s(s + a) 

[e + Sgn e] 

L 

T^ is a transcendental function of 
a 2 A/MK c (Fig. A-5 and Ref. A-5) 

K c /s 2 

[e + / 2MK c |e| sgn e] 

K c s 3 

{e + (l/3)e 3 + Wee + w[(l/2)e 2 + We] 3/2 } 
W = +1 for [e + (l/2)e|e| ] > 0 
= -1 for [e + (l/2)e|e| ] < 0 
( Ref . A-6 ) 
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Figure A-5. 

Lead Compensation at Switching Point 
For Controlled Element Y = K /s(s + a) 
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The decision logic model for the time-optimal control algorithm be- 
haves like a function switch and accounts for the initial increase in the 
time delay (beyond that due to quasilinear tracking) in response to a step 
input • 

The model of Fig, A-4 should thus serve as one possible explanation of 
operator behavior in response to random plus step inputs. 
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APPENDIX B 


EFFECTS OF HORIZONTAL HIND SHEAR GRADIENT 
ON HEIGHT CONTROL WHILE TRAVELING 
IN "DOLPHIN" MANEUVERS 


The effects on closed-loop characteristics of a headwind-to-tailwind 
shear gradient, 3u g / 3h < 0, with decreasing altitude are shown in 

Figs. B-l through B— 4. The airspeed regulation subsidence, 

wJJ = 0.3 rad/sec, corresponds with a velocity-to-pitch attitude control 
high-frequency gain K u = -0.01 rad 9 c /ft/sec in both sets of figures. 
Figures B-l and B-2 represent the effects with the lower value of height— 
to-collective control high-frequency loop gain = 0.0245 in/ft/sec 
(damping ratio 0.785; undamped natural frequency 0.51 rad/sec); Figs. B-3 
and B-4 represent the effects for a higher value of = 0.0923 in/ft/sec 
(damping ratio 0.7; undamped natural frequency 1.0 rad/sec). The shear 
gradient 3u g /3h varies along the complex and real root loci in Figs. B-l 
and B-3 as shown in the corresponding logarithmic Sigma-Bode root loci. 

Figs. B— 2 and B— 4. The results in Figs. B— 1 through B— 4 are expanded 

further in this appendix to represent the effects of two other values of 
airspeed regulation subsidence, oj” = 0.15 and 0.6 rad/sec, corresponding 
respectively to two other values of velocity-to-pitch attitude control 
high-frequency gain, viz., K u = -0.005 and -0.02 rad 9 c /ft/sec. 

Figures B-5 through B-8 present the results for 0 )^ = 0.15 rad/sec, and 
Figs. B-9 through B-12 present the results for = 0.6 rad/sec. 

Not surprisingly, the higher height control bandwidth in Figs. B-3 and 
B-4 (u£ = 0.3 rad/sec), Figs. B-7 and B-8 ( to” = 0.15 rad/sec), and 

Figs. B— 1 1 and B-12 (coJJ = 0.6 rad/sec) provides a higher threshold of 
instability caused by the wind shear gradient. Although the height con- 
trol bandwidth in Figs. B-3, B-4, B-7, B-8, B-ll, and B-12 is only doubled 
with respect to that in Figs. B-l and B-2 (a/J = 0.3 rad/sec), Figs. B-5 
and B-6 (u£ = 0.15 rad/sec), and Figs. B-9 and B-10 (u)[J = 0.6 rad/sec), 

the wind shear gradient required to cause instability is between 2-1/2 and 

3 times greater in Figs. B-3, B-4, B-7, B-8, B-ll, and B-12 depending on 
the speed regulation subsidence. 
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joj (rad/sec) 



Real s 


UH-1H at 20 kts, Sea Level, 8000 lbs, raid c.g. , = 0.0923 in/ft/sec 

-Z w = 0.567 rad/sec -Zg * 9.29 ft/sec 2 -in -X = 0.005 rad/sec 

c u 

" z u = °* 198 rad/sec (X a -g)K u = u£ = 0.15 rad/sec 1/T. = 1,13 rad / sec 

K u = -0.005 rad 9 c /ft/sec 


Figure B-7 . Complex Root Locus of Characteristic Values for the Closed-Loop 
Determinant of Height and Airspeed Control Equations While Traveling Showing 
Effects of Headwind-to-Tailwind Shear with Decreasing Altitude 9u / 9h < 0 


Complex Conjugate Roota^frcm Height Regulation 
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Figure B-8. Sigma-Bode Root Locus Showing Effect of Headwind-to-Tailwlnd Shear 

with Decreasing Altitude: 3u /3h < 0 




UH-1H at 20 kts, Sea Level, 8000 lbs, mid c.g., = 0.0245 in/ft/sec 


-Z w = 0.567 rad/ sec 
-Z u = 0.198 rad/sec 


-Z* = 9.29 ft/sec^-in 

c 

(X -g)K =co" =0.6 rad/sec 
u 

K u = -0.02 rad 9 c /ft/sec 


-X u = 0.005 rad/sec 

1/Tr = 1.13 rad/sec 
L h 


Figure B-9. Complex Root Locus of Characteristic Values for the Closed-Loop 

Determinant of Height and Airspeed Control Equations While Traveling Showing 

Effects of Headwind-to-Tailwind Shear with Decreasing Altitude 3u /3h < 0 
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Complex Conjugate Roots from Height Regulation 
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Figure B-12. Sigma-Bode Root Locus Showing Effect of Headwind-to-Tailwind Shea 

with Decreasing Altitude: 3u,,/3h < 0 



If the speed subsidence is smaller (u)” 0.15 rad/sec and 0.3 with 
u)[j* 1.0 rad/sec) than the height control bandwidth, an increasing wind 
shear gradient of the sense presented causes a monotonic reduction in 
height control bandwidth but a slight increase in height control damping 
ratio before it causes loss of damping and instability. This is a more 
favorable circumstance than if the speed subsidence approaches (o^ =0.3 
with ■ 0.5 rad/sec and w” - 0.6 with wfj * 1.0 rad/sec) or exceeds 

(a)JJ = 0.6 with = 0.5 rad/sec) the height control bandwidth, whence an 
increasing wind shear of the sense presented causes simultaneous and mono- 
tonic loss of height control bandwidth and damping ratio. In all cases 
shown, the wind shear causes a monotonic increase in the speed subsidence. 

In Figs. B-5 and B-6 (wJJ = 0.15 rad/sec) (lower closed-loop height 
control bandwidth and lowest speed subsidence), a shear gradient of 
-0.2/sec (-11.84 kt/100 ft) decreases the closed-loop height bandwidth 

about 36 percent and increases the closed— loop damping ratio about 

12 percent. Figures B-5 and B-6 also show that the shear gradient must 

increase to -1.155/sec (-68.39 kt/100 ft) to cause instability. 

The larger closed-loop height bandwidth shown in Figs. B-7 and B-8 
U” = 0.15 rad/sec) (coupled with the lowest speed subsidence) provides 
the greatest margin of stability against this particular type and sense of 
shear gradient among the six examples shown. A shear gradient of 

-1.31 /sec (-77.57 kt/100 ft) reduces the closed-loop height bandwidth 
about 44 percent and increases the closed-loop damping ratio about 

21 percent. Figures B-7 and B-8 show that the shear gradient must 
increase to -3.78/sec (-224 kt/100 ft) to produce instability. 

Consider next the two examples with the intermediate values of the 
speed subsidence o>” ■ 0.3 rad/sec. In Figs. B-l and B-2 

( a)” = 0.3 rad/sec) (lower closed-loop height control bandwidth), a shear 
gradient of only -.0842/sec (-5 kts/100 ft) decreases the closed-loop 

height bandwidth about 20 percent without altering the damping ratio. 
Figures B-l and B-2 also show that there is still a considerable margin of 
stability vested in the damping ratio, because a shear gradient of 
-1.125/sec (-66.6 kts/100 ft) is required to produce instability. 
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The larger closed-loop height bandwidth shown in Figs. B-3 and B-4 
(w(| = 0.3 rad/sec) provides a greater margin of stability against the 
headwind-to-t ail wind shear with decreasing altitude. A shear gradient 
3Ug/3h of -0.85/sec (-50 kts/100 ft) reduces the closed-loop height 
bandwidth by about 33 percent without altering the closed-loop damping 
ratio significantly; whereas, a shear gradient of -3.22/sec 

(—190 kts/100 ft) is required to produce instability. Consider finally 
the two examples with the largest value of the speed subsidence, 
“u = 0*6 rad/sec. 

In Figs. B-8 and B-9 (oj^ = 0.6 rad/sec) (lower closed-loop height 
control bandwidth), the shear gradient causes simultaneous and monotonic 
decrease in closed-loop height bandwidth and damping ratio. For example, 
a gradient of -0.193/sec (-11.43 kt/100 ft) reduces the height bandwidth 
about 19 percent and reduces the damping ratio about 24 percent; whereas, 
a shear gradient of -1.09/sec (-64.5 kt/100 ft) is required to produce in- 
stability. This example provides the smallest margin of stability against 
this particular type and sense of shear gradient among the six examples 
shown. 

The larger closed-loop height bandwidth shown in Figs. B-ll and B-12 
K = 0.6 rad/sec) provides a greater margin of stability against the 
shear gradient, although the decreases in height bandwidth and damping 
ratio are again simultaneous and monotonic with increasing gradient. A 
shear gradient of -0.568/sec (-33.51 kt/100 ft) reduces the closed-loop 
height bandwidth about 27 percent and reduces the damping ratio about 
14 percent; whereas, a shear gradient of -2.64/sec (-156 kt/100 ft) is 
required to produce instability. 

The results from these six examples are summarized in Figs. 26 through 
28 in the main body of the text. 
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APPENDIX C 


HEIGHT CONTROL IN THE BOB-DP OR -DOWN RAPID RESPONSE PHASE 


Height Control With Vertical Velocity 
Command-Height Hold (VCHH) 


During the rapid response phase of the vertical maneuver with velocity 

command-height hold capability in Fig. 32 (p. 75) in the main body of the 

text, the reference for the compensatory height regulation loop via the 

controller is slewed by the velocity command integration while the height 

feedback loop remains closed to provide a measure of immunity against 

disturbances— especially the destabilizing effect of a headwind-to- 

tailwind shear gradient with decreasing altitude: 3u /3h < 0. 

§ 

Concurrently, the three other compensatory loops for longitudinal (Fig. 32 
in the text) and lateral stationkeeping (Fig. 33, p. 76 in the text) and 
heading-hold are closed via the controller throughout the rapid response 
phase of the vertical maneuver. To simplify the subsequent analysis, we 
shall consider only the vertical and longitudinal motions by assuming that 
the controller decouples the lateral-directional motions during the 
vertical maneuver by means of the methods described, for example, in 
Ref. C-l . 

The closed-loop transfer function for height response, h, to a verti- 
cal velocity command, h c , can be expressed with good approximation at low 
frequencies which are typical of the guidance bandwidth as in Eq. C-l. 


s) 


h 

c 



Closely Coupled Stationkeeping Dipole 


C-l 



where 



<* * t ~) - z j K i;( 

C 



3u 

'6 ‘'h ;o “6 "V‘L “ Th* Z u 

c c h 


+ iT-- 


( C— 2 ) 


4 l^“xl 4 s2 + (V« V 

T Lh A TJ in Fig. 32 

T t A t. in Fig. 32 

l, x x 



Not of much practical significance 

because (X u - /Zg Z ) is so small 
u c u c 

for rotorcraft 


K* = K t T* in Fig. 32; K. = K r T. in Fig. 32; 
h x x 


K. A 

= K z in Fig. 32 

L x 

x 


K -\K lnFls - 

c he 



K* A 

K in Fig. 32 

L h h 


to provide a unit gain for the closed-loop velocity response to a velocity 
command at low frequencies in Eq. C-l in the absence of wind shear gra- 
dient 3u /3h. 

o 



to provide matched lead and lag frequencies in the absence of wind shear 

gradient 3u /3h. Therefore, 

6 




c 


1 


Z 6\ 

c h 


if 3u / 3h = 0. 
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The corresponding closed-loop transfer function for the collective 
control displacement, 5 co n » required by a vertical velocity command, h c , 
can be expressed with good approximation for the purpose of guidance as in 


Eq. C-4 


'"'coll, „ 

--- - - (s) 
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N? C ° n 
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A" 


where 




K. (s + 
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tT 

L h 

c 



8u 
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9h 




Closely coupled 

Stationkeeping 

Dipole 



(C-4) 


Example Topic (a) Expanded from Main Text 
Beginning on Page 74 


Example of design for a special case of VCHH with preservation of 
inherent heave damping , i.e., 1 /k: A i/t * !/ T = 1/T fl and no wind 

shear, i.e., 3u g /3h = 0. Therefore, -Z fi K£ = -Zg = 1/T Q , where 

1/Tq^ is the inherent heave damping, a function of main rotor properties, 
disc loading, geometry and true airspeed, and Z { , the collective control 
effectiveness, Is also a function of main rotor properties, disc loading, 
geometry, and true airspeed* The closed— loop height response, h, to 
vertical velocity command, h Q , will be from Eq. C-l : 



(C-5) 


closely coupled 
dipole 


closely coupled 
stationkeeping 
dipole 
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The closed-loop collective control required will be from Eq. C-4: 


N ?coll 

S coll r . A h c . . . 

■ (s) = — (s) = - 


s ( s + + F“)Ux ;a) x^ 

1 _ . e 2 9 2 

W s ( s , + tt )(s * 
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closely closely closely coupled 
coupled coupled stationkeeping 

dipole dipole dipole 


Equation 05 approximates the controlled element form Y c = K c /s(s + a) 
with K c = a for which the generic time-optimal velocity command and rapid 
height response to a step command in height has already been described in 
Figs. 34 and 36 the main text (p. 78 and 80, respectively). Equation 06 
shows that the required collective control displacement will be 
proportional to and in phase with the vertical velocity command time 
history as described in Fig. 35 in the main text (p. 79). 


Example Topic (e) Expanded From Main Text 
Beginning on Page 88 


Example of redesign of VCHH system to increase tolerance to a wind 

shear gradient 8u„/ 9h - -1.0/sec (-59.2 kt/100 ft). Redesign and T t 
S h 

in Eq, 02 so that 

1/T *2 ' * 2/T »2 a ” < ’ ' 3 V 3h Z “ * < l / T 9 2 )2> 

if 3Ug/3h = -1.0/sec 

Z u = -0.198/sec 

-Zr = 9.29 f t/sec 2 -in 

°c 

and 1 /Tq = 0.567/sec 

0 2 

therefore, -Z, K: * 1/T fl and 1 /T t = 0.198 T ft + 1/T fi = 0.916/sec. 
o c n 0 2 L h a 2 9 2 

If -Zj ± 1 /Tl = 1/Tg as in Example Topic (a), the transfer 
c c h 2 

function for height response to vertical velocity command (Eq. C-l) in the 
presence of the designated wind shear gradient will now be 
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(C-9) 
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A N h 
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T I s + T" J 
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; ( S+ T^“) ( S . + T^~) 



closely coupled 
dipole In wind 
shear 


closely coupled 
stationkeeping 
dipoles 


Equation C-9 Is the same as the height response in Eq. C-5 and Fig. 36 in 
the text (p. 80) for Example Topic (a), except that the vertical velocity 
command in Fig. 34 in the text (p. 78) is optimized for the effect of this 
particular wind shear. 


The transfer function for closed-loop collective control required in 
the presence of the designated wind shear gradient will now be 


M ?coll 

S „ A N h 
coll, x c , N 

- (s) = — (s) = 

’ A" 



dipole in wind 
shear 


(C-10) 


stationkeeping 

dipole 


where -1/Zf T Q 
^c 0 2 


0.0610 in/ft/sec. 


Equation C-10 is, however, different from the collective response required 
in Eq. C-6 and Fig. 35 in the text (p. 79) for Example Case (a), because 
C— 10 represents the cost of controlling collective displacement to 
provide a time optimal height response in the designated wind shear with 
VCHH. Figure 40 (p. 85) shows the collective time history required by 
Eq. C-10. 
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Example Topic (f) Expanded from Main Text 
Beginning on Page 88 


Redesigned Example Case (e) in the absence of wind shear 
(3Ug/8h * 0). Since the redesign in Example Case (e) is based on such an 

extreme value of the wind shear gradient, it is necessary also to examine 
the height and collective response to a vertical velocity command in the 
absence of wind shear. If 9u / 3h = 0 with the VCHH design gains and 

5 

lead/lag compensation selected in Example Case (e), Eq. C-9 no longer 
represents the height response; instead, Eq. C-2 for changes to 

have the following values 


% 


1 

-Z. K h 

c n _ n „ 2(0.567) 

0.721/sec and ? h 2(0.721 > 


9u 

0.787 if — ■ & = 0 
9n 


and Eq. C-l , in turn, changes to have the values in Eq. C-ll for the 
height response to a vertical velocity command 


s) 


h 

c 



0.567 (s + 0.567) U';uf] 


s[0 .787 ;0 .721 ] 


/ 

closely coupled 
stationkeeping 
dipole 


(C-ll) 


The effective controlled element in Eq. C-ll is still approximately of the 

. 

form Y c « K c /s(s + a), where a = 0.862/sec based on the dimensionless 

"k • 

variability and a = 1.0/sec based on the half-power frequency; but the 


* M Bandwidth" is a vague term unless the displayed signal spectrum is 
rectangular. For other spectral shapes, the dimensionless variability can 
be used to define a rectangular bandwidth equivalent, i.e., 


LC * 11 <“> d “] 2 

<D. = — £ — X (Ref * C ~ 2) 

e L da) 
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steady-state gain of the velocity response to a velocity command is less 
than unity; (it is 0.618). The corresponding collective response transfer 
function in Eq. C-12 is based on Eq. C-4 with appropriate modification of 
[? h ;u) h ] " [0 . 787 ;0 . 72 1 ] . 


coll 


5 ii 

„.coll 
4 N h 


(s) = 


A" 


<s) - - 


S(S + V )(s 

9 2 0 2 


Z_ T 


6 c 9 2 s [0.787 ;0 .721 ] 

X X 

/ 


closely coupled 
stationkeeping 
dipole 


(C-12) 


where -1/Zg T 0 = 0.0610 in/ft/sec 
1 /Tq^ = 1/Tq^ = 0.567/sec 

If the high frequency gain (-Z^ K£ = 0.567/sec) is increased by a 
factor of 1/0.618 = 1.618 to compensate for the reduced gain in Eq. C-ll, 
the new high frequency gain in Eq. C-lla is 0.917, and the corresponding 
high frequency gain in Eq. C-12a is 0.0987 


1 .618N* 

Ns, h = 


A" 


0.917(s + 0.567) ] 

(s) = — — 2 — £_ 

s[0.787;0.721 ] [c";oo”]- 

X X 


(C-lla) 


and 


coll 


Closely Coupled Stationkeeping Dipole 


<s) = 



(C-12 a) 
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If the amplitude and time intervals of the velocity command are recomputed 

to account for the increased high frequency gain, K , and increased 

2 

effective heave damping, a, in the maneuver urgency factor K c M/a A = 1.0, 
where M = |h c |, then 

|h | - a 2 A/K c = (0.862) 2 50/0.917 = 40.51 ft/sec, 

where a = 0.862 based on the dimensionless variability 

T g - 1.61/1.0 = 1.61 sec 
T f = 0.511/1.0 = 0.511 sec 
and T c = 2.12 sec 

where a = 1.0/sec based on the half-power frequency of the vertical 

# • 

velocity (h) response to a velocity command (h c ). The application of the 
two different values for a is deliberate because of the form of the 
controlled element in Eq. C-lla. 

The modified velocity command is shown as a function of time in the 
main text. Fig. 43 (p. 89). Transient responses of collective and height 
displacements for this modified (and nearly re-optimized) velocity command 
are given in the maintext Figs. 44 and 45 (pp. 90 and 91), respectively. 
Corresponding velocity response is presented in the phase plane, Fig. 46 
(p. 92), which confirms that the height and vertical velocity responses 
have been nearly re-optimized in the absence of wind shear. 

Conclusion from these examples with VCHH 

The "worst-case" design for the wind shear in Example Case (e) is 
preferable, because it results in a conservative design which can be more 
easily re-optimized as in Example Case (f) by velocity command amplitude 
and time interval adjustment in the absence of wind shear. 
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APPENDIX D 


THE LONGITUDINAL SPEED RESPONSE TO PITCHING MOMENT CONTROL 
OR PITCH ATTITUDE CONTROL FOR ARTICULATED 
ROTARY-WING AND TILT-ROTOR AIRCRAFT 


In this appendix, we shall derive the closed-loop longitudinal speed 
response for typical rotorcraft based on approximations for the ratio of 
moduli and the ratio of real parts of the complex dipole in the longitu- 
dinal speed response to pitching moment control for articulated rotary- 
wing and tilt-rotor aircraft. The derivation will be performed for 

hovering flight, because the simplicity of hovering dynamics will enable 
us to illustrate all of the essential features of the complex dipole which 
represents the stabilized "pendulum effect" without introducing the un- 
necessary algebraic complications caused by higher forward speed. 
Nevertheless, the results are valid at forward speeds typical of formation 
station-keeping • 

When aerodynamic stability derivatives which are usually small are 

omitted, the two longitudinal equations of hovering motion which describe 
the pendulum effect in attitude and speed can be written (Refs. D-l 

through D-4) together with a third pitch attitude command and feedback 
equation with rate augmentation for stabilization as: 


s - X 

u 
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"V 
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-X 
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-M s(s 
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0 

► = { 

-M 

. u + 
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0 K s 
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CD 
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6 b 


0 

j g 

*e) 


0 

c 


where 

g = acceleration of gravity 


Kq = pitching rate feedback gain to the longitudinal 
pitching moment control in dimensional units of 
time 



K0 * dimensionless pitch attitude feedback gain to the 

longitudinal pitching moment control 

M * pitching moment applied to the rotorcraft 

M q = (l/I y )(3M/3q) 

\ = (l/ly)(3M/3u) 

- (l/I y )(3M/35 B ) 

s = complex Laplacian operator 

X = longitudinal force applied to the rotorcraft 

X u » (l/m)(3X/3u) 

X 5 - (l/m)(3X/3 5 ) 

8 B 

I y = pitching moment of inertia of the rotorcraft 
m ■ mass of the aircraft 

q * pitching rate 

u * longitudinal perturbed velocity 

Ug = longitudinal atmopsheric gust velocity 

= pitching moment control displacement 
0 = perturbed pitch attitude 

3 = partial differential operator 


The characteristic equation of this set is called the (stabilized) "hover- 
ing cubic" and is given by 


s(s - H )(s - X ) + gM 


(y + y [Mj (» - Xj + X s Mj 


V ^ / 

"Hovering cubic" 
of vehicle alone 


Contribution from Pitch 
Attitude Stabilization 
with Rate Augmentation 
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Although the speed stability derivatives Xy and My can be rather 
large, they generally have only a small effect on the stabilized pendulum 
effect in attitude control in the absence of external disturbances such as 
longitudinal gust velocity, Ug. Consequently we can further simplify the 
stabilized hovering cubic by neglecting the small products involving the 
speed stability derivatives and factoring the characteristic equation 
approximately into a quadratic representing the stabilized pendulum effect 
and a subsidence representing a reduced surge. damping, 0 < 1/Tg < (-Xy) 


4 ' . [s 2 - M s + M { K . + M { Kg] (s + 1/Tg ) 

M B H B 1 

■ [s 2 + K K „ - Mjs + Mj Kg] (s + 1/Tg ) 

B n ^ B 1 

= [s 2 + 2? p u p s + (o) p ] 2 ] (s + 1/T'^) 


where twice the real part of the complex pole is 2? u = Mr K - M and 

, * P P °B ^ <1 


/ n * 

the square of its modulus is (u> 1 - M* K rt . 

P « B 9 


For articulated rotary wing and tilt rotor aircraft, the corresponding 
speed response numerator to the pilot's <5g input in the absence of an 
external disturbance, u , is:* 

O 



Xx (s^ + 2? ci) s + 
o v u u 


*The corresponding speed response numerator to the pilot's pitch 
attitude command input, 9 C , is 

"9,. : VSg (s 2 + 2t u » u s + u 2 ) 

D-J 



where twice the real part of the complex zero is : 


25 <*> 
u u 


-M 


and the square of its modulus is : 




_ JL 

X r 


M. 


At low forward speeds and in hovering, u>- 



therefore, 



K 0 


expresses the ratio of squared moduli of the complex dipole representing 
the stabilized pendulum effect. Typically -Kq is on the order of unity. 
The ratio of real parts of the complex dipole is 


5 w 
u u 

t t 

5 

P P 


-M 


M, K - M 

6 b q q 


Since < 0 and > l M ql> the ratio of real parts will usually be 

less than unity; and the amplitude of the residual oscillatory pendulum 
component in longitudinal speed perturbations will be noticeable (and may 
be bothersome) to the pilot when he regulates disturbances via 0 C through 



We can represent the effect of speed regulation by the pilot with the 
aid of the block diagram below 


u 

g 



The transfer function for the attitude-controlled rotorcraft with pitch 
stability and control augmentation system (SCAS) has been derived here as 

K 0 X 6 [s 2 + 2C to s + w 2 ] 

B _ u u u _ 

(s + 1/T' ) [s 2 + 25'u's + U') 2 ] 

9 1 P P v P J J 

where K 0 - -1, = g » 32.2 ft/sec 2 , and 0 < 1/Tg^< -Xy 



EXAMPLE 

An effective rotor response delay of 0.1 sec will be added to a 
pilot's effective delay of 0.33 sec so that t u = 0.43 sec by assumption. 
Typical values for the other singularities are listed below. 


0.0433 
2.1 rad/sec 

/ 

Up =2.48 rad/ sec 


<1/T 0 )' - 0.0257 rad/sec c u = 
5p - 0-337 u u = 
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Thus the open-loop transfer function for speed regulation by the pilot 
will be 



-32.2K u [0.0433;2.1]e" 0,43s 
(0.0257) f 0.337 ;2.48] 


where the following abbreviated notation is used for polynomial factors in 
root locus form: 


Real factor (■=!— ) means (s + 



Quadratic factor [?;w] means [s^ + 24 ms + 

A Bode diagram of this open-loop transfer function is shown in 
Fig. D-la attached. Notice the decade-and-a-half (from 0.03 to 
1.0 rad/sec) of K/s controlled element characteristics which are favorable 
to unit gain crossover with ample phase margin. A corresponding complex 
root locus is shown in Fig. D-lb attached. Closed-loop speed regulation 
characteristics are identified for a typical range of the pilot's gain 
0.0223 < -i^ < 0.06 rad/ft/sec in speed regulation tasks demanding high 
precision. 

u 

The rotorcraft s airspeed response numerator, Nx , provides favorable 

B 

second-order lead compensation at 2.1 rad/sec, which is nearly equal to 
the stabilized "pendulum” frequency at 2.48 rad/sec resulting from the 
pitch attitude closure with rate augmentation in the SCAS. Consequently 
the pilot's pure gain airspeed loop closure in Figs. D-la and D-lbaccom- 
plishes two desirable results: (1) it increases the closed-loop 
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Figure D-1Tj. Complex Root Locus for Airspeed. Outer-Loop Closure, 

u — ►» 3_ with Pitch AFC 3 * * 




"pendulum" bandwidth in excess of 3 rad/sec while maintaining an adequate 
damping ratio in excess of 0.4, and (2) it restores an oscillatory closed- 
loop "phugoid" with a relatively high undamped natural frequency, 
1.5 rad/sec, which is relatively invariant over a wide range of damping 
ratio, from 0.9 to 0.3, -as the pilot's gain is increased. The closed-loop 
"phugoid" governs the airspeed regulation bandwidth, but will appear to 
the pilot using a precision display of speed as an oscillatory mode with 
increasing amplitude as he increases his gain. The airpseed regulation 
bandwidth is not sensitive to a two-fold variation in the pilot's gain in 
excess of -Ky * 0.0223 rad/ft/sec, which variation primarily affects damp- 
ing ratio and modal response coefficient. For example, if we hypothesize 
a step change in commanded airspeed, the transient overshoot associated 
with the closed-loop "phugoid" will be nil if the pilot's gain 
-Ky = 0.0223 rad/ft/sec, but the overshoot will increase from 6 percent to 
28 percent as the pilot's gain, -Ky, is increased from 0.03 to 
0.06 rad/ft/sec. Thus the net effect of the complex dipole in the speed 

response is to cause a marked loss in precision of airspeed regulation as 

. 

the pilot's gain is increased above = 0.03 rad/ft/sec, even though the 
bandwidth is relatively invariant. 

Up to this point in the example we have illustrated the effects of a 
pure gain strategy for airspeed regulation by the pilot. We can next 
represent the pilot's trimming strategy for speed control by an integral 
gain, Kjj/s, in parallel with the pure gain, Ky, both of which operate on 
speed error, u e , as shown in the following block diagram. 





The open-loop transfer function for speed regulation will now become, 
again using the abbreviated notation for factors, 



K 


W« 

D U 




where K^qX^ - -32. 2K^ as before 

From this transfer function it is clear that successful trimming by the 
pilot requires that - 1/T^. Physically this means that the pilot's 

trimming gain ratio, 1^/1^, must match the reduced surge damping, 
0 < 1/T^ < (-Xy). In terms of the Bode diagram of Fig. D-1^ this will 
extend the K/s controlled element characteristic to zero (to provide zero 
steady-state speed error) and create a real closed-loop dipole in the 
vicinity of I^/l^ - 1/T^, which represents the trimming subsidence. 
Thus, except for this dipole near the origin and the pole at the origin, 
the complex root locus will be virtually identical to that in Fig. D-lh. 

This example can be extended again to illustrate the effects of a 

hovering strategy by the pilot which is intended to regulate longitudinal 

position, x. A block diagram for this part of the hovering task is shown 
below. 
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The open-loop transfer function for longitudinal position regulation will 
become, again using the abbreviated notation for factors, 



K 

K K fl X«. C=£)[5 jm 
u 0 o t , v K ' l u u- 
B u 


-T s 
u 



where K^jKqX^ 


B 


-32.21^ as before 


Note the similarity in form to the previous open-loop transfer function 
for trimming speed regulation. In the present instance, however, the 
pilot's strategy for his displaceraent-to— rate gain ratio, operating 
on position error will depend on his relative thresholds of indifference 
to velocity and displacement cues and his perception of the disturbance 
environment. Three types of position regulation strategy are depicted in 
the sketches of complex root loci in the following figures. 

In Fig. D-2 , the disturbance environment. is presumed to be relatively 
benign, K^/Ky > 1 /Tq^, and, except near the origin, the complex root locus 
is virtually indistinguishable from that in Fig. D-lb. 

In Fig. D-3 the disturbance environment is presumed to be signifi- 
cant, Kjj/Ky » 1 /Tq^ , and the predominant closed-loop "phugoid" bandwidth 
remains proportional to the pilot's gain, -K^, with a relatively re- 
stricted variation in closed-loop damping ratio about the optimum value, 
0.707. These closed-loop characteristics enable the pilot to suppress the 
disturbances as long as his achievable bandwidth (approximately 1 rad/sec) 
exceeds that of the disturbances. Otherwise his gain (and bandwidth) will 
necessarily regress. 

In Fig. the pilot is presumed to have adopted too high a gain 
ratio, Kjj/Kjj, which results in a lower than desirable closed-loop damping 
ratio for the "phugoid," regardless of his gain, -K^. 
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Figure D-3 • Complex Root Locus for Longitudinal Position Loop Closure 
x, u 9 in a Disturbed Environment with Gain Ratio K JYL 
Nearly Optimum (k/k^=0.3) a U 
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APPENDIX E 


NUMERICAL SOLUTIONS OF THE TRANSCENDENTAL INVARIANCE 
CONDITIONS FOR TIME OPTIMALITY 


Table E — 1 provides numerical solutions of the transcendental invar- 
iance conditions for time optimality given in Table A-l in Appendix A for 
controlled element Y £ = K c /s(s + a). Definitions of symbols are as 
follows : 


2 

a A/K c M Controlled element-and-maneuver characteristic (dimensionless 

inverse maneuver urgency factor) 

a Characteristic damping, rad/sec 

A Amplitude of step function representing command input, units of 

displacement appropriate to the maneuver 

K c High frequency gain of the controlled element in appropriate 

units 


M 


Average absolute amplitude of the controller response appro- 
priate for each element assuming the controller response to be 
bang-bang with equal positive and negative amplitudes 

Starting pulse time interval for the controller response, sec 

Final pulse time interval for the controller response, sec 

Time to complete the rapid response maneuver, T g + T^, sec 

Characteristic lead compensation time at the pulse switching 
point of the controller, sec 
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TABLE E-l. INVARIANCE CONDITIONS FOR TIME 
OPTIMALITY WITH A STEP INPUT AND CONTROLLED ELEMENT 

Y c • K c /s(s + a) 


a 2 A/K c M 


aT, 


T f - aT p 


T s + T f 
• T c - aT c 


T, 


aT T 


0 

• 

0 

• 

• 

0 

• 
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0 

• 

• 

CO 

9 

• 

• 

• 

• 

✓e“ . 

fT 

o 
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2 3 

T T 

S , S 

T s " 2! + 3! 

e o 

2x f 

/J 

3T f 

2 3 

T T 

S , s 

E o " 2! + 3! 

• 

• 

• 


• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

• 

0.01 

0.118913 

0.05769 


0.17660 

35.015 

0.03 

0.21056 

0.09975 


0.31031 

20.395 

0.10 

0.40194 

0.18091 


0.58285 

11.398 

0.25 

0.674 

0.281 


0.955 

7.39 

0.50 

1.03 

0.387 


1.42 

5.53 

0.75 

1.33 

0.458 


1.79 

4.73 

1.0 

1.61 

0.511 


2.12 

4.27 
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0.583 


2.73 
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• 

• 

• 

• 

• 

+ T f 

• 
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• 

• 

e + 2Zn2 

OO 

• 
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(Continued) 
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TABLE E-l 


(Concluded) 



e 


o 





fr 

o 




0.01 

1.19 

0.03 

1.22 

0.1 

1.27 

0.25 

1.35 

0.50 

1.46 

0.75 

1.54 

1.0 
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1.5 
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• 
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• 

• 

• 

e 

• 

e + t 

oo f 

00 

/e 


CD 
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0.476 
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£ n 2 



00 


1.77 
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1.65 
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3.70 

1.77 

2.01 
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4.10 
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2.12 
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2.26 

2.23 
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2.57 
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• 

• 
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• 

• 
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• 

/e 

oo 

m 

e + 2 in 2 
00 


(1 - £ n 2 ) 

4(1 - JLn 2 ) 


00 
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Controlled Element -and -Maneuver Characteristic (dmls) 

a 2 A/K c M 

(More Aggressive-^ (inverse Maneuver Urgency Factor) ►-(Less Aggressive) 

Lead Compensation at Switching Point 
For Controlled Element Y = K /s(s + a) 

C C' ' 
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APPENDIX F 


EQUATIONS AND METHODS USED FOR VARIANCE ANALYSIS 

A. LONGITUDINAL AND VERTICAL 

A block diagram of the mathematical models used for variance estima- 
tion in height control and longitudinal control tasks is given in 
Fig* F— 1 . The pilot's control strategies and divided attention noise 
sources are represented at the left side of Fig. F-l . Command equaliza- 
tion for velocity control with position hold capabilities are represented 
across the upper and lower portions of Fig. F-l. Atmospheric turbulence 
noise sources are shown right of center, and rotorcraft dynamics and kine- 
matics are depicted at the right side of Fig. F-l . 

Laplace-transformed linear differential equations of motion for the 
longitudinal and vertical variance analysis are given in Table F-l, to- 
gether with symbols not otherwise identified in the block diagram. 
Fig. F-l. 

B. LATERAL AND DIRECTIONAL 

A block diagram of the mathematical models used for variance estima- 
tion in directional control and lateral control tasks is given in 
Fig. F-2. The pilot's control strategies and divided attention noise 
sources are represented at the left side of Fig. F-2. Control equaliza- 
tion for yaw rate command with heading hold capability is represented 
across the upper portion of Fig. F-2, and control equalization for lateral 
velocity command with position hold capability is represented across the 
lower portion of Fig. F-2. Atmospheric turbulence noise sources are shown 
right of center, and rotorcraft dynamics and kinematics are depicted at 
the right side of Fig. F-2. 

Laplace-transformed linear differential equations of motion for the 
lateral and directional variance analysis in hover are given in Table F-2, 
together with symbols not otherwise identified in the block diagram. 
Fig. F-2. 
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Figure F-1. Block Diagram of Mathematical Models for Covariance Estimation 
in Height Control and Longitudinal Control Tasks 







TABLE F-l. LAPLACE-TRANSFORMED LINEAR DIFFERENTIAL EQUATIONS OF MOTION FOR 
LONGITUDINAL AND VERTICAL COVARIANCE ANALYSIS 
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TABLE F-l (CONCLUDED) 


SPECIAL SYMBOLS 


[N] 

2 6 

" s ~ T s 

*4 

where 

[D] 

2 6 

= s + — s 

T 

X 

+4 

where 


TRIMMED STABILITY AND CONTROL DERIVATIVES 
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Figure F-2. Block Diagram of Mathematical Models for Covariance Estimation 
in Directional Control and Lateral Control Tasks 




















TABLE F-2. LAPLACE-TRANSFORMED LINEAR DIFFERENTIAL EQUATIONS OF MOTION FOR 
LATERAL AND DIRECTIONAL COVARIANCE ANALYSIS 
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TABLE F-2 (CONCLUDED) 


SPECIAL SYMBOLS 

[N] = s 2 - s + where x = x in Line 5 

x z v 

[D] - s 2 + s + where x = x in Line 5 

x z ^ 

TRIMMED STABILITY AND CONTROL DERIVATIVES 


Y 

v 



N 

v 





or x 


y 


in Line 
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or x in Line 6 
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